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HIGH-PRESSURE FLAME SYSTEM FOR POLLUTION STUDIES WITH 


RESULTS FOR METHANE-AIR DIFFUSION FLAMES 

Irvin M. Miller and Howard G. Maahs 
Langley Research Center 


SUMMARY 

A high-pressure flame system has been designed and constructed for studying 
nitrogen oxide formation in fuel-air combustion. Its advantages and limitations 
have been demonstrated by tests with a confined laminar methane-air diffusion 
flame over the pressure range from 1 to 50 atm. The methane issued from a 
3.06 -mm-diameter port concentrically into a stream of air contained within a 
20. 5 -nim- diameter chimney. Except for special studies of flame stability, fuel 
and air flow rates were held constant at Ul.8 standard cm^/min (seem) and 
2450 seem, respectively. 

As the pressure of the flame was increased from 1 to 50 atm, its shape 
changed from wide and convex to slender and concave. At the same time, there 
was a marked increase in the amount of carbon in the flame as evidenced by its 
greatly increased luminosity. The height of the flame changed moderately with 
pressure. At higher pressures the limits of stability and extinction were 
greatly reduced and the limits within which the flame smoked were greatly 
increased. 

Average nitric oxide formation rates were greater than those predicted by 
the conventional two-step Zeldovich mechanism with oxygen atoms assumed to be 
in equilibrium with oxygen molecules . Without making the oxygen-atom equilib- 
rium assumption, required superequilibrium oxygen-atom ratios were calculated 
using the Zeldovich mechanism. At 1 atm, these ratios were close to literature 
values reported for premixed flames, but they increased markedly above 5 atm. 

It may be that at elevated pressures an increasing amount of nitric oxide was 
formed through nitrogen-hydrocarbon reactions in addition to the Zeldovich 
mechanism. 

The molar emission index for N0 X (the moles of total nitrogen oxides 
formed per mole of methane consumed), which was O .85 x 10“3 a t 1 atm, increased 
rapidly as pressure was increased; it reached a maximum value of 1.21 x 10 “^ at 
about 9 atm; and thereafter it decreased slowly until at 50 atm it was back to 
the same value that it was at 1 atm. This maximum occurred in spite of the 
fact that the peak temperature in the flame decreased from about 2130 K at 
1 atm to about 1850 K at 20 atm, and thereafter, at higher pressures, decreased 
relatively little. Unfortunately, the present data are insufficient to permit 
a clear choice among several alternate explanations of this maximum offered in 
this report. A number of experiments which would aid in making such a choice 
are suggested. 



INTRODUCTION 


Pollution from aircraft turbine engines is an increasing national concern 
as more aircraft and more powerful engines are placed into service. The prob- 
lem is particularly severe in and around airports because of aircraft concen- 
tration and the greater quantity of fuel burned on take-off. Of the various 
pollutants discharged from these engines, two of the most harmful are nitric 
oxide and nitrogen dioxide, commonly referred to as NOy. A primary national 
goal is to control these NOy emissions, and numerous studies, both experimental 
and analytical, have been made toward this end. 

Experimental studies have involved measurements of pollutants in the com- 
bustors and exhausts of turbine engines to correlate engine emissions with 
parametric adjustments in the combustor. Analytical combustion models have 
been formulated to predict pollutant formation in a given combustor. Each 
model consists of a unique combination of assumed chemical and fluid-mechanical 
mechanisms which is supposed to represent the combustion process. Such models 
can be of great value for designing more efficient, cleaner combustors if the 
assumed combination of mechanisms accurately represents the combustion, flow, 
and pollutant formation processes actually occurring in a given combustor. 

Since the actual combustion and flow processes are very complex (involving fuel 
injection, atomization and evaporation, air addition, partial recirculation, 
diffusion, heat transfer, turbulence, etc., at pressures up to 30 atm), experi- 
mental verification is necessary. However, such verification is difficult and 
very expensive in actual combustors. An alternate approach would be to check 
combustion models with a simplified laboratory combustion system. To the 
authors’ knowledge, there are no such laboratory combustion systems being used 
to verify combustion models at high pressures (30 at m and higher). 

Within an actual turbine engine combustor, the mixing state of fuel and 
air can exist locally between two limits: completely unmixed, as in the vicin- 

ity of a vaporized fuel droplet near the fuel jet, and completely mixed, as in 
the highly turbulent region of the combustor. These two localized limiting 
states can be simulated in the laboratory by either of two simplified combus- 
tion systems, the data from which can be useful in verifying current combustion 
models as well as in furthering the understanding of the NOy formation process 
in hydrocarbon flames. These two systems are a laminar diffusion flame and a 
laminar premixed flame. The laminar diffusion flame was chosen for the initial 
design and construction of the experimental facility because it is experimen- 
tally safer and simpler, and it is not complicated by problems of mixing and 
flashback which become severe for premixed flames at high pressures. 

On the basis of the foregoing considerations, a high-pressure flame system 
was designed, constructed, and tested to determine the effect of pressure on 
NOx emissions from hydro carbon- air flames. Descriptions of this facility, of 
the combustion gas analysis system, and of the instrumentation for measuring 
the temperature of the luminous zone in the flame are presented in the first 
three sections of this report. Results obtained for a methane-air diffusion 
flame are presented and discussed in the last two sections. 
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Certain commercial materials are identified in this report in order to 
specify adequately which materials were used in the research effort. In no 
case does such identification imply recommendation or endorsement of the pro- 
duct by NASA, nor does it imply that the materials are necessarily the only 
ones or the best ones available for the purpose. In many cases equivalent 
materials are available and would probably produce equivalent results. 
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SYMBOLS 

molar concentration of species i, moles/cm3 
average molar concentration of species i, moles/cm^ 
mass concentration of carbon, g/cm^ 
first Planck constant, 3*74l844 pj-cm 2 /sec 
second Planck constant, 1.438833 cm-K 

average carbon-particle diameter, A or ran (1 A = 0.1 nm) 

Planck function, defined by equation (D7) 
function defined by equation (D3) 

emission index, (g of N0 X as if all were N0 2 )/(kg of CH^ consumed) 

molar emission index, ^moles of N 0 X produced)/(mole of CH 4 consumed) 

molar emission index for NOx at equilibrium, 

(moles of N 0 X produced)/(mole of CHj^ consumed) 

equilibrium constant defined by equation ( 10 ) 

spectral extinction coefficient, cm 2 /g 

arbitrary constant 

forward and reverse elementary reaction-rate constants, cm^/mole-sec 
(number in subscript (e.g., kpg) refers to equation number) 

radiation path length, cm 

molecular weight of species i, g/mole 

complex index of refraction, n - inK 

Avogadro's number, 6.023 * 10 2 ^ molecules/mole 

number density of carbon particles, cm” 3 
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number of carbon atoms available in methane fuel per unit volume of 
reaction zone in flame assuming stoichiometric combustion in flame 
zone 

refractive index 

molar flow rate, moles /sec 

pressure, atm (l atm = 101.3 kPa) 

universal gas constant , cal/mole-K or cm^-atm/mole-K 
reaction rate of species i, moles/cm^-sec 
average reaction rate of species i, moles /cm^-sec 
stoichiometric molar fuel -air ratio 
temperature, K 

adiabatic flame temperature, K 
average temperature in reaction zone, K 
brightness temperature at wavelength A, K 
peak temperature in flame, K 
time, sec 

o 

volumetric flow rate, cirr/sec 
reaction volume, cm^ 

specific surface reaction rate, moles /cm^-sec 
mole fraction of species i 
average mole fraction of species i 

mole fraction of methane that would exist in flame zone under 
stoichiometric conditions if none had been converted to soot 

equilibrium oxygen-atom mole fraction for stoichiometric premixed 
combustion 

oxygen- atom mole fraction required by equation (8) to give observed 
NO^ formation rate 

dispersion exponent; arbitrary constant 
arbitrary constant 


h 



e A 

K 

X 

p c 

<J> 

4 > f 
x 


surface area of carbon particles in flame per unit volume of reaction 
zone in flame 

spectral emissivity 

absorption index 

wavelength, A or ni (1 A = 0.1 nm) 
carbon-particle mass density, g/cm^ 
equivalence ratio 

pseudoequivalence ratio, defined by equation (l) 
quantity defined in appendix E 


Subscripts : 

max maximum value assuming all carbon in fuel has been converted to 

carbon particles 

stoich stoichiometric conditions 

STP at standard temperature and pressure, 293 K and 1 atm 

Superscripts : 

carbon for stoichiometric conditions assuming that all methane has been 

transformed to carbon and hydrogen 

E experimental 

stoich for stoichiometric conditions assuming perfect mixing 

Chemical symbols: 


C 

c 2 

CH 

CH 4 

CO 

co 2 

Ho 


h 2 o 


carbon atom 
diatomic carbon 
carbon-hydrogen radical 
methane 

carbon monoxide 
carbon dioxide 
hydrogen molecule 
water 
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I 


He 

helium 


M 

arbitrary third body 


N 

nitrogen atom 


n 2 

nitrogen molecule 


NO 

nitric oxide 


no 2 

nitrogen dioxide 


N0 X 

total nitrogen oxides 

(NO 

0 

oxygen atom 


°2 

oxygen molecule 


OH 

hydroxyl radical 


Abbreviations : 


C/L-A 

chemiluminescent analyzer 

cal. 

calibration 


diam 

diameter 


F.M. 

flowmeter 


IR 

infrared 


i . d. 

inside diameter 


o.d. 

outside diameter 


seem 

standard cm3/min 



DESCRIPTION OF SYSTEM 

The high-pressure flame system developed for this study represents the 
first such system for the investigation of confined laminar diffusion flames at 
elevated pressures up to 50 atm. The confined laminar diffusion flame, first 
analyzed by Burke and Schumann (ref. l) , is a flame burning at the exit of a 
tube from which fuel gas issues coaxially into a cylindrical duct of coflowing 
air (or other oxidizer). An important part of the present high-pressure system 
is the provision for collection and analysis of the products of combustion from 
the flame, particularly oxides of nitrogen. The initial design of the system 
was based on the following parameters: 
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Fuel Methane 

Molar fuel-air ratio 0.1 to 1.0 

Fuel flow rate, seem 152 maximum 

Combustion pressure, atm 1 to 50 


Burner Design 

The heart of the high-pressure diffusion flame system is the burner, where 
methane gas (99*97 percent pure) issues coaxially from a cylindrical tube into 
a concentric cylindrical quartz tube through which air (99*995 percent pure, 
total hydrocarbon concentration less than 1 ppm) is coflowing, as shown in fig- 
ure 1. If the air and fuel streams are laminar, then the air and fuel mix 
largely by diffusion. Under these conditions, a stable diffusion flame results 
upon ignition. Burke and Schumann (ref. l) were the first to analyze theoret- 
ically this type of flame, and they developed an equation which predicts the 
flame height and shape for laminar flow. However, they experimentally verified 
their equation only over the very limited pressure range from 1 to 1.5 atm. 

The equation developed by Burke and Schumann (ref. l) and Reynolds criteria 
for laminar flow in tubes were used to size the burner for a small luminous 
flame. In addition to laminar flow requirements, the flame had to be luminous 
and nonsmoking. A luminous flame was necessary for measuring the temperature 
of the luminous zone by a relatively simple spectroscopic technique. This tech- 
nique is discussed in a subsequent section. The requirement that the flame be 
nonsmoking avoids the problem of contaminating the system as well as of removing 
soot from the combustion products. 

To obtain a smoke-free, luminous flame at 1 atm, the pseudoequivalence 
ratio cf) ! , defined as the actual molar ratio of fuel to air divided by the 
stoichiometric molar ratio of fuel to air,^- must be less than 1. If <J> 1 is 
near 1, the flame is luminous and sooty, and if <J) 1 is greater than 1, the 
flame is nonluminous and contains large quantities of unburned hydrocarbons. 
Examples of these types of flames are shown in figure 2. (Extensive, detailed 
descriptions of such flames for the butane-air system at 1 atm are given by 
Barr (ref. 2).) From appendix A, the value of <J> 1 can be set by controlling 
the volumetric flow rates of fuel and air at standard conditions of temperature 
and pressure (293 K and 1 atm), or 


1 v STP,fuel 
S V STP,air 


( 1 ) 


The term ”pseudoequi valence ratio” has been coined for the present dif- 
fusion flame because the concept of equivalence ratio does not hold the same 
fundamental significance for a diffusion flame that it does when used for pre- 
mixed flames. In an overventilated diffusion flame (i.e., (j) 1 < l) such as the 

one under study, it is generally assumed that the effective equivalence ratio 
in the flame zone where the reaction actually occurs is very close to unity. 
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where S is the stoichiometric fuel-air ratio and Vgrpp f ue p and Vgipp a -j_ r 
are the volumetric flow rates of fuel and air at standard conditions. For all 
measurements, the volumetric flow rates of fuel and air were maintained at 1*1.8 
and 2 1*50 seem, respectively, for reasons discussed in the section. Smoke, 
Stability, and Extinction Limits. Under these conditions with S = 0.105 for 
methane, <j> ? = 0.162 by equation (l). 

The fuel tube has an inside diameter of 3.06 mm and a wall thickness of 
0.05 mm. (See fig. 1.) The outlet edge of the tube is tapered to minimize the 
formation of turbulent eddies. The quartz air tube, hereafter referred to as 
the quartz chimney, has an inside diameter of 20.5 mm. In order to obtain a 
flat velocity profile in the fuel and air streams as they exit from the fuel 
and air ducts (with the intent of simplifying any theoretical analysis of the 
combustion process in the diffusion flame), a porous disk of sintered stainless 
steel was installed near the exit of the fuel duct as well as near the exit of 
the air duct, as shown in figure 1. In preliminary tests, no porous disks were 
installed in the fuel and air ducts, and it was not possible to sustain a flame 
after ignition at pressures approaching 50 atm. At these pressures, with the 
hot-wire igniter located near the fuel duct exit, the flame appeared and dis- 
appeared in a cyclic manner - on for a fraction of a second, off for more than 
a second. However, with porous disks installed in both the fuel and air ducts, 
a flame, once ignited, was self-sustaining. Other tests with a porous disk 
only in the fuel duct were partially successful - combustion was sometimes 
cyclic, other times continuous. Consequently, all subsequent tests employed 
porous disks in both the fuel and air ducts. 

In early tests at pressures above 40 atm, a band of tiny condensate drop- 
lets formed on the inside of the chimney, extending from about halfway up the 
flame to the top of the chimney. This was probably due to the cooling effect 
of the nitrogen pressurizing gas sweeping past the side of the chimney. This 
problem was considerably reduced by surrounding the chimney with a larger square 
quartz chimney. As a result, only a very thin band of condensate formed inside 
the chimney above the flame, and this soon reached a steady-state amount. 


Sample-Handling System 

When methane is burned with air, the primary combustion products are 
carbon dioxide and water vapor, although small quantities of carbon monoxide 
and NOx are also formed. To make studies of pollution formation in the flames, 
a sampling system was devised which collects the total combustion products 
given off by the flame and transports them to the gas analyzers. A description 
of the gas analyzers is given in a subsequent section. 

A schematic diagram of the sample collection system is shown in figure 3. 
The hot combustion products from the flame exit from the top of the quartz 
chimney and are swept upward into an electrically heated total-sample collector 
by the nitrogen pressurizing gas. The diluted combustion gases continue to 
flow upward through an electrically heated tube, a hot-air-heated 5-ym particle 
filter, and a hot-air-heated back-pressure regulator. The filter and regulator 
are located outside the pressure chamber. At the exit of the back-pressure 
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regulator, the total gas sample is vented to unheated sample transfer lines at 
atmospheric pressure and thence to the gas analyzers. All surfaces contacted 
by the total gas sample, from the chimney exit to the gas analyzers, are of 
stainless steel. 

The purpose of the nitrogen dilution gas is twofold: (l) it serves as 

the pressurizing gas in the combustion chamber, and ( 2 ) it reduces the partial 
pressure of water vapor in the total sample and hence lowers the saturation 
temperature of the water vapor in that sample. This makes it possible to main- 
tain all the water produced during combustion in the vapor state by heating the 
sample transfer lines to only 393 K, even for the highest pressure of 50 atm. 
For the sample transfer lines downstream of the back-pressure regulator where 
the sample is at 1 atm, no heating is necessary. (See appendix B.) 


Ignition System 

Ignition is accomplished by means of an electrical resistance-heated 
hot-wire igniter. This igniter, shown in figure 4, is retractable through a 
Chevron packing gland at the top of the sample transfer line in the chamber 
(see fig. 3). During operation, the igniter is positioned 10 to 20 mm above 
the burner, and electric power is applied until ignition is attained. Power to 
the igniter is then shut off and the igniter retracted to the top of the sample 
transfer line where it is completely out of the sample flow. 


Flow and Pressure Control System 

A schematic diagram of the flow and pressure control system is shown in 
figure 5* Gas flows are controlled with micrometer needle valves. The upstream 
pressure of these valves is controlled by dome-loaded regulators and the down- 
stream pressure or chamber pressure is controlled by a dome-loaded back-pressure 
regulator. The control pressures for the three dome-loaded regulators and for 
the back-pressure regulator are provided by two separate spring-loaded regula- 
tors. Since the spring-loaded regulators tend to drift, block valves are used 
to lock in the set pressures. A gas reservoir was installed in each locked-in 
system to minimize the effect of leaks and changes in ambient temperature. 

Pressure in the pressure chamber is measured by two separate Bourdon tube 
gages. Low pressures from 1 to 8 atm are measured with a gage having a pressure 
range from 0 to 1.03 MPa (150 psia) and an accuracy of 0.066 percent of full 
scale. High pressures from 8 to 50 atm are measured with a gage having a pres- 
sure range from 0 to 6.89 MPa (1000 psig) and an accuracy of 0.1 percent of full 
scale. Flow rates of methane, air, and nitrogen are measured with linear ther- 
mal mass flowmeters having an accuracy of 1 percent of full scale. Meter ranges 
were: 0 to 152 seem for methane, 0 to 3000 seem for air, and 0 to 50000 seem 

for nitrogen. The methane and air flowmeters were calibrated with a soap-film 
meter whose calibration is traceable to the National Bureau of Standards. Dur- 
ing calibration, appropriate corrections were made for barometric pressure and 
ambient temperature. Calibrations were made at 1, 25, and 50 atm, the meters 
being zeroed at each pressure prior to calibration. 
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Sight Glass Windows 


Safety sight glass windows of 66.7-mm viewing diameter are located on 
opposite sides of the pressure chamber and are made to withstand a pressure of 
50 atm. The burner and chimney can be viewed through these windows. Each 
flange-mounted sight glass consists of two 25-mm-thick quartz plates laminated 
with a sheet of Teflon (FEP), type A, nominally 0.0381 mm thick, compressed 
between the plates to a thickness of 0.0127 mm. 

Color photographs of the flame were taken through these windows with a 
35-mm single-lens reflex camera fitted with a 200 mm/f2.8 reflective lens and 
with color negative film. Proper exposure of the flame at each pressure was 
determined by trial and error. 


Operation 

The high-pressure flame system is operated by setting the desired flows 
and pressure in the pressure chamber and igniting the flame. The pressure 
upstream of the micrometer needle flow control valves is first set at a suit- 
able pressure above the desired ambient pressure in the chamber and the flow- 
meters are zeroed. Then the desired chamber pressure is set with the back- 
pressure regulator and the chamber is pressurized with nitrogen (99-995 percent 
pure). When the desired chamber pressure is reached, fuel, air, and nitrogen 
flows are adjusted to their target values: Ul.8 seem for methane, 2^50 seem 

for air, and about 6000 seem for nitrogen. (The nitrogen flow rate of 6000 seem 
is high enough to prevent condensation in the sample transfer lines but not so 
high as to dilute excessively the combustion gases to be analyzed. No flame 
instability was created by the nitrogen flow. ) The electrical heaters for the 
sample transfer lines are energized, and after the sample lines reach the 
desired temperature, the igniter rod is lowered to a point close to the burner, 
the flame is ignited, and the igniter rod is withdrawn from the stream of com- 
bustion products. 


CHEMICAL ANALYSIS 

All chemical analyses are made on the total gas mixture exiting from the 
high-pressure chamber: combustion products, excess air, and pressurizing 

nitrogen. This overall, or total, gas analysis approach was selected instead 
of a detailed profile analysis approach for several reasons. First, local com- 
position measurements by probe sampling would have presented severe experimental 
difficulties since the sample probe and its associated manipulating mechanism 
must be located entirely within the enclosed high-pressure chamber and actuated 
remotely. Secondly, reaction zones typically become thinner as pressure is 
increased, so that there was some question as to whether satisfactory spatially 
resolved measurements could be obtained in any event. Finally, for these 
studies, it was believed to be adequate to measure the average pollutant con- 
centration generated by the flame as opposed to measuring detailed concentra- 
tion profiles throughout the flame. 
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In early checkout tests with the high-pressure system, all products were 
analyzed by gas chromatography with the exception of NOx, which is not amenable 
to gas chromatographic determination at the low concentrations encountered. 

The analysis for NOx was performed with a commercial chemiluminescent NO-NOx 
analyzer having a lower limit of detectability of about 0.3 ppb . Chemical 
analyses from these checkout tests for overventilated flames (the proposed 
normal operating condition) showed CO to be below the limit of detectability 
(about 100 ppm) and the total concentration of hydrocarbons was typically about 
6 ppm, indicating that combustion was virtually complete to CO2 and H2O. When 
the system was operated at near-stoichiometric flow rates of CH^ and air, the 
expected tall smoky flame was encountered. This condition was avoided so as 
not to deposit soot within the chamber or sample lines. However, several tests 
were intentionally made on underventilated flames. These showed substantial 
quantities of CO as well as CO2, and also substantial quantities of unburned 
fuel along with four additional hydrocarbons detected by flame ionization. 

These hydrocarbons were tentatively identified as acetylene, ethylene, ethane, 
and methanol; however, no attempts were made to make a more positive identifi- 
cation because underventilated flames were not of concern in this study, and in 
any case, are customarily avoided in most practical combustion systems. 

Since these early tests showed that for the proposed normal overventilated 
operating condition, virtually all the fuel carbon appeared as CO2, consider- 
able simplification of the analytical system could be made by substituting a 
continuous infrared CO2 monitor for the gas chromatograph. Analysis by infra- 
red absorption is not only easier than gas chromatography ,- but also a continu- 
ous record is obtained. The concentration of CO2 is of particular interest 
because it provides a means of determining the extent of nitrogen dilution in 
the pressure chamber independent of measurements of the separate fuel, air, and 
nitrogen flow rates (these measurements are inherently less accurate). For all 
subsequent experimental tests, the concentration of CO2 was determined with a 
(commercial) infrared analyzer having a repeatability of ±1 percent. 

For the measurement of NO and NOx concentrations, a (commercial) chemi- 
luminescent N 0 -N 0 x analyzer was employed. Corrections for quenching by CO2 
and H2O in the sample were made according to the factors reported by Maahs 
(ref. 3 ). Because of substantial dilution of the combustion products by the 
pressurizing N2 , these corrections were fairly small, about 2.5 percent. Con- 
version efficiency of NO2 to NO in the NOx catalytic converter (employing a 
stainless steel coil) was measured frequently and was typically above 98 per- 
cent. Appropriate corrections were made for the actual conversion efficiency 
and are reflected in the total NOx concentrations. 

A schematic diagram of the final gas analysis system is shown in figure 6. 
The bypass flow around the chemiluminescent analyzer was a convenient means of 
providing adequate sample flow to the analyzer while at the same time insuring 
the sample pressure to the analyzer to be 1 atm. This precaution was taken 
since the chemiluminescent analyzer is mildly sensitive to sample pressure. 

The flowmeter at the exit of the infrared CO2 analyzer was installed to monitor 
the flow rate through the infrared cell, which in effect, is an extremely sen- 
sitive way to monitor cell pressure. By maintaining this flow rate constant, 
constant cell pressure is insured. The other vents and bypass flow lines were 
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installed as a convenient means of insuring that switching between sample and 
calibration gases did not disrupt flows in the overall analytical system. 


Details of the gas chromatographic system employed in the early checkout 
tests are given in appendix C. Although sample analyses by gas chromatography 
were made only in these early tests, the column arrangement and operating con- 
ditions might be of some interest to the reader since this gas chromatographic 
system proved to be not only convenient to use but also permitted an analysis 
for water and at the same time completely avoided the need for backflushing any 
of the columns. 


FLAME TEMPERATURE MEASUREMENT 

The measurement of flame temperature is by far the most difficult, most 
uncertain, and least satisfactory of any of the measurements made. Measurement 
of temperature in a flame is inherently difficult, but the poor physical acces- 
sibility of the present flame for probing with thermocouples is the main problem. 
Consequently, optical methods had to be used. Sodium line reversal was ruled 
out for several reasons. Not only is the technique complicated, but also, no 
convenient way could be devised to introduce the sodium into the flame without 
encountering flameout and/or sodium deposition throughout the chamber and 
sample lines. Furthermore, it was felt best to avoid the introduction of for- 
eign materials (i.e., sodium) if at all possible so as not to perturb the NOy 
formation rates. Refractive index methods such as those described by Weinberg 
(ref. 4) were tried, but distortions caused by the double chimneys and the 
5. 0-cm-thick quartz windows were too great to permit meaningful analysis. 

Considerable effort was expended in attempting to apply the self-absorption 
method described by Kuhn and Tankin (ref. 5 ), hut without success. In this 
method, optical techniques are used for determining the absorption coefficients 
of the carbon particles naturally occurring in the flame. Kuhn and Tankin 
demonstrated the feasibility of the method for a large (12.5 cm high) luminous 
propane-air diffusion flame issuing from a 2 . 5-cm-diameter burner. But for the 
present, less luminous, much smaller methane flame (about 10 mm high, issuing 
from a 3 . 06 -mm-diameter burner), the necessary sensitivity could not be obtained 
to get a reasonable estimate of the carbon-particle absorption coefficients. 
Another closely related method also attempted is that of Kurlbaum as described 
by Ribaud, Laure, and Gaudry (ref. 6). It was not successful either, largely 
for the same reasons. 

Only the exceedingly tedious method of OH absorption offered any hope of 
being able to measure accurately local flame temperatures, and even this method 
is somewhat uncertain because of the small thickness of the combustion zone in 
the flame and the consequent high resolution required of an optical system. 
Therefore, it was decided that a measurement of a spatial average temperature 
near the hottest point in the flame was the best that was reasonably possible. 

For this measurement , the method of two-color pyrometry was employed with the 
naturally occurring carbon in the flame as the radiant emitter. In applying 
this method, it is assumed that the spectral emissivity of the carbon par- 
ticles can be expressed in the form 
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where k is an arbitrary constant, L is the path length, X is the wave- 
length, and a is the dispersion exponent. Since the exponent in the exponen- 
tial term is typically much smaller than unity for most flames, the exponential 
term was approximated by the first two terms in a Taylor series expansion, giving 
ct 

z kL/X . Then, by using Planck’s black-body law and the definition of emis- 
sivity, and by making brightness temperature measurements T-^ ^ and T *^ \2 
at two wavelengths X^_ and X 2 , k can be eliminated and the flame tempera- 
ture T^ can be expressed by 
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where c 2 is the second Planck constant, 1.438833 cm-K. The only nonmeasured 
parameter required to calculate flame temperature from this equation is the 
dispersion exponent a, taken for the present methane-air flames to be 1.4. 
McAdams (ref. 7) recommends an average value for a of 1.39 for all hydro- 
carbon flames, while d'Alessio and coworkers (ref. 8) report a value of 1.4 for 
methane-oxygen flames. However, T f is actually relatively insensitive to the 
precise value selected for a. An error analysis of equation (3) using reason- 
able values for X^ and T-^ ^ indicates that a 30-percent uncertainty in a 

produces only about a 3-percent uncertainty in T^. 

Brightness temperature measurements were made with two disappearing- 
filament optical pyrometers fitted with interference filters at different wave- 
lengths , nominally 5000 A and 6500 A. Since the wavelengths of these filters 
are known with great accuracy, ^ the largest error in the flame temperatures 
calculated by equation (3) is due to experimental errors in measuring T-^^. 

To minimize these errors, as well as the error propagation in equation (3) 
caused by the wavelengths being too close together, initial tests were made 
with the pyrometers sighting on a flame through a number of different optical 


^The 5000-A filter has a central wavelength of 5002 A, 63-percent trans- 
mission, and 83-A bandwidth at half-height. The 6500-A filter has a central 
wavelength of 6508 A, 60-percent transmission, and 90-A bandwidth at half- 
height . 


13 


filters at wavelengths ranging from 4000 A to 7000 A. Ultimately, 5000-A and 

o 

6500-A filters were selected as the best compromise between wavelength separa- 
tion and eye sensitivity. Both pyrometers were calibrated for brightness tem- 
perature as follows: The true temperature of a temperature-controlled black- 

body cavity was established by a standard optical pyrometer calibrated by the 
National Bureau of Standards. The pyrometers to be calibrated were then 
sighted on the black body through the same laminated quartz sight glass window 
as was used in the chamber and through two thicknesses of quartz glass repre- 
senting the walls of the two quartz chimneys surrounding the burner. Readings 
were made on the standard and test pyrometers over a suitable temperature range 
For making actual experimental measurements, the pyrometers were attached to 
adjustable mounts having both vertical and horizontal motion capability and 
were sighted through the windows in the pressure chamber on the flame at the 
position of interest. Typically, this position was just below the tip of the 
flame in the region of the most intense radiation. The experimental pyrometer 
readings were then converted to standard black-body brightness temperatures by 
the calibration curves, and flame temperature was calculated from equation (3). 

During the course of the experiments, the pyrometer system was improved 
by incorporating both filters into a single pyrometer with provisions for 
switching between the two filters. In this way, possible uncertainties in 
positioning two individual pyrometers when viewing the flame were effectively 
eliminated. 


RESULTS AND DISCUSSION 

The high-pressure flame system was successfully operated over the pressure 
range from 1 to 50 atm. All tests were made with methane as the fuel and air 
as the oxidizer. Results from these preliminary tests revealed a number of 
interesting and sometimes unexpected phenomena regarding flame shape, carbon 
formation, smoke emission, flame stability, extinction limits, NOx emission 
levels, and flame temperature. These phenomena are discussed in the following 
sections . 


Flame Shape and Carbon Formation 

A very marked change occurred in the shape and structure of the methane- 
air diffusion flame as pressure was increased from 1 to 50 atm, as illustrated 
in the color photographs in figure 7. All photographs were taken at the same 
flow rate of methane (1+1.8 seem) and air (2450 seem), and are shown at an 
enlargement of about 4.6 to 1. The bottom edge of each photograph in figure 7 
is coincident with the top edge of the burner, and the lateral position of the 
burner is shown below each photograph. At 1 atm, the sides of the flame bow 
outward from the inner fuel core into the surrounding air annulus, with the 
result that the diameter of the flame is considerably larger than the diam- 
eter of the fuel tube. The blue zone in the photographs is likely due to O 2 
Schumann-Runge radiation and CO + 0 continuum radiation (see Gaydon and Wolfhard 
(ref. 9)), and the brighter blue-green or green zones on the sides and near the 
base of the flame are likely due to C 2 and CH radiation (ref. 9). A zone of 
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•weak carbon-particle luminosity is evidenced by the off-white zone in the 
upper portion of the flame. As pressure was increased above 1 atm, the out- 
ward bowing of the sides of the flame continually decreased. The most nota- 
ble change in shape occurred below about 5 atm. At about 10 atm the flame 
exhibited a slight concavity and the diameter decreased considerably. Also, 
the zones of molecular radiation became narrower and did not extend as high; 
the carbon zone became yellower and much brighter. A further decrease in flame 
diameter continued as pressure was increased above 10 atm; however, most of the 
change in shape and structure occurred when 20 atm was reached. At higher pres- 
sures up to 50 atm, only minor changes occurred. It seems likely that a change 
in reaction mechanism with increasing pressure, as evidenced by the presence of 
increasing amounts of carbon (a reaction intermediate) in the flame, is largely 
responsible for the pronounced change in shape of the flame. Because of their 
relatively high mass, solid carbon particles, once formed, do not diffuse out- 
ward from the gaseous fuel core as does the gaseous fuel itself, and combustion 
must proceed by virtue of the fact that oxygen diffuses inward to the carbon 
zone. Hence, the flame assumes its constricted appearance. 

Interestingly, the height of the flame did not change drastically, at 
least up to about 30 atm, although a rather marked change occurred in the shape 
of the flame with pressure. A plot of flame height as a function of pressure 
is shown in figure 8, where the flame heights were determined from the color 
photographs in figure 7. The height to the tip of the observable blue reaction 
zone and the height to the tip of the luminous carbon zone are separately indi- 
cated for pressures at which they differ. On the basis of the position of the 
luminous tip, there is roughly a 25 -percent increase in the flame height between 
1 and 10 atm, and on the basis of overall height, a decrease of about 10 per- 
cent. Above 10 atm, the flame height decreases slowly up to about 30 atm; then 
it drops more steeply to its lowest value at 50 atm. With the exception of 
50-atm flames, the flame height was typically 10 mm. 

These flame height data in figure 8 conflict with theoretical predictions 
that the height of a confined laminar diffusion flame should be independent of 
pressure at constant mass flow both for equal fuel and air velocities (see 
refs. 1 and 10) and for unequal fuel and air velocities (see refs. 11 and 12). 
Experimental verification of these theoretical predictions has been previously 
demonstrated only at very low pressures from 0.06 atm to 1.5 atm. In contrast, 
Parker and Wolfhard (ref. 13) experimentally showed for acetylene that inde- 
pendence of flame height with pressure at constant mass flow is true only for 
very small flames; for larger flames the height is independent of pressure only 
in those situations where diffusion is not the sole rate-limiting process. 

Larger flames were stated to increase in height as pressure increases, particu- 
larly when carbon is formed in the flame. Differing somewhat from these obser- 
vations in reference 13 are the theoretical predictions of Powell (ref. l4) for 
a "sandwich burner," which consisted of a row of alternating fuel and air slots. 
For constant mass flow and for pressures at which diffusion is controlling, 
Powell predicts flame height to be independent of pressure; but at very low 
pressures at which reaction rate is controlling, he predicts flame height to 
vary inversely with pressure. None of these predictions adequately describe 
the observed change in flame height of the present methane-air flame over the 
entire pressure range from 1 to 50 atm. But none of the theoretical predictions 
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account for carbon formation in the flame either, and it is believed that this 
increasing amount of carbon formed with increasing pressure is largely respon- 
sible for the experimental disagreement with theoretical predictions. 

In addition to their prediction of the independence of flame height with 
pressure, Burke and Schumann (ref. l) similarly predict that the shape of the 
flame should be independent of pressure. This is obviously not indicated by 
the present high-pressure results, nor even by low-pressure results. (See 
refs. 11, 12, and 13.) However, an accurate prediction of flame shape may be 
too much to expect from the admittedly simple theory of Burke and Schumann 
which employs many assumptions including those of constant properties (imply- 
ing, among other things, no buoyancy forces), infinite reaction rate, and no 
carbon formation; that the theory predicts flame heights as well as it does for 
small changes in pressure is indeed remarkable. 

The observed increase in carbon formation with increasing pressure in 
a diffusion flame is a well-known phenomenon. (See the literature review of 
McArragher and Tan (ref. 15).) At 1 atm, very little carbon is formed as indi- 
cated by the very low carbon luminosity of the flame. (See fig. 7.) As pres- 
sure increases, the amount of carbon formed increases considerably to produce a 
brilliant incandescent flame. The brightness of these higher pressure flames 
cannot be appreciated from the photographs in figure 7, since camera exposure 
time had to be reduced considerably from 3 sec at 1 atm to 0.0005 sec at 50 atm 
so as not to overexpose the film. Some idea of this increase in luminosity o can 
be seen from figure 9 where measurements of brightness temperature at 6500 A 
made in the tip of the flame have been converted to monochromatic emissive power 
at 6500 A by Planck 1 s law and are plotted against pressure. The fuel and air 
flow rates were maintained at the low constant levels mentioned previously to 
avoid unoxidized carbon (or soot) from being liberated from the flame as smoke 
at any of the test pressures from 1 to 50 atm. 

Figure 7 shows that carbon formation (as evidenced by the white or yellow 
luminous zone) occurs mainly in the tip of the flame at low pressures, and as 
pressure is increased, the carbon zone moves downward to fill an increasingly 
larger portion of the flame. A rough estimate of the actual mass concentra- 
tion c of carbon in the tip of the flame can be obtained from experimental 
data on brightness temperature and a knowledge of the flame temperature itself 
(presented subsequently in this report). Details of the calculation are given 
in appendix D. The mass concentration of carbon in the flame is shown in fig- 
ure 10 as a function of pressure. Carbon concentration increases roughly 
3 orders of magnitude from 2.5 atm to 20 atm, and above that, not much further 
increase occurs. In appendix D, the relationship between the number density of 
carbon particles Nq and the mass concentration c is also derived. If rea- 
sonable values of average particle diameter dQ are assumed to be on the order 
of 50 A to 100 A ( 5 run to 10 nm) (see Narasimhan and Foster (ref. l6)), then 
curves of particle number density as a function of pressure can be constructed, 
as shown in figure 11. To obtain a rough estimate of particle number density 

o 

for a 1-atm flame, the particle number density curve for 50-A-diameter particles 
in figure 11 was extrapolated from 2.5 atm down to 1 atm to what appears to be 
a reasonable, though admittedly arbitrary, value of 10® particles/cm®. Fig- 
ure 11 illustrates that by far most of the increase in carbon-particle number 
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density occurs between 1 and 5 atm, which is in agreement with what would be 
surmised intuitively from the photographs in figure 7. Above 5 atm, there is 
more than an order of magnitude increase in particle number density up to about 
20 atm, above which the number density changes relatively little on up to a 
pressure of 50 atm. This phenomenon of increasing carbon concentration with 
increasing pressure is of particular concern because it clearly indicates that 
the mechanism of hydrocarbon oxidation at high pressures in diffusion flames is 
markedly different from that at low pressures. Hence any mathematical descrip- 
tion of the diffusional combustion of methane (or any other hydrocarbon, for 
that matter) with air at pressures even moderately above 1 atm, must properly 
account for the formation and subsequent combustion of carbon particles in the 
flame . 


Smoke, Stability, and Extinction Limits 

Similar to carbon formation, smoke formation also increased with increasing 
pressure. The terms "carbon" and "smoke" are used here to distinguish between 
carbon particles within the flame (as evidenced by flame luminosity) and uncom- 
busted carbon particles liberated from the flame. It was possible at 1 atm to 
burn methane at a fairly high flow rate, producing a large smokeless incandes- 
cent flame, but a small increase in pressure transformed this smokeless flame 
into a quite smoky one. Also, at higher pressures, the limits of stable combus- 
tion were considerably reduced. To define more clearly these regions of smoke 
and stable combustion, a limited series of tests were made at 1, 5, 20, and 
50 atm for the present methane-air diffusion flame over a wide range of fuel and 
air flow rates. The results are shown in figure 12. The main purpose of these 
tests was to determine regions of extinction and stable combustion, with a 
secondary purpose of locating luminous and nonluminous regions and smoke bound- 
aries. The upper limits of 152 seem for methane flow rate and 3000 seem for 
air flow rate in figure 12 were dictated by the upper limits of the flowmeters 
used. No attempt was made in these studies to produce an extensively detailed 
map of flame shape and extinction such as that produced by Barr (ref. 2) for a 
butane-air flame; the boundaries indicated in figure 12 are approximate only. 

A number of pertinent observations can be made from figure 12. Perhaps 
the most striking of these is that the limits of combustion are considerably 
reduced as pressure is increased; in other words, at higher pressures the flow 
rate ranges of methane and air over which combustion can be maintained are much 
smaller than those at lower pressures. It was also noted that extinction of a 
flame burning at high pressures was often abrupt, whereas at lower pressures 
considerable flickering and instability occurred in a flame burning near its 
extinction limit. This is believed to be related to the fact that reaction 
rate increases as the pressure increases, while diffusion rate remains rela- 
tively constant with pressure; consequently, the reaction zone becomes thinner. 
The higher pressure flames, with their characteristically thin reaction zones, 
are apparently more susceptible to a sudden extinction by a slight perturbation 
because the reactants cannot be transported into the thin reaction zone fast 
enough to reestablish a flame which is burning under marginal conditions. At 
lower pressures, slight perturbations on the flame tend to create only flicker- 
ing and instability because the transport of reactants is sufficiently fast 
relative to the reaction rate to bring reactants into proximity of the thicker 
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reaction zone. Sufficient combustion is thus maintained somewhere within the 
zone, and the flame reignites itself throughout the reaction zone in hursts. 

Figure 12 also shows that the region of stable combustion is markedly 
reduced as pressure is increased. For the present purposes, stable combustion 
was taken to be indicated by a virtually motionless flame without any flicker- 
ing, wavering, pulsating, or other transient phenomena. At 1 atm the flame is 
everywhere stable within the combustion boundaries shown, except at low fuel 
flow rates near the transition boundary between an overventilated and under- 
ventilated flame. At 5 atm (fig. 12(b)) the region of stable combustion is 
already reduced considerably in size and occurs only at relatively small fuel 
flow rates. At 50 atm the stable region is very small indeed (fig. 12(d)), and 
regardless of flow rate, all flames at 50 atm were found to be luminous. Even 
at only 20 atm, a nonluminous flame could be produced only in a very limited 
region of low fuel and air flow rates. 

On the basis of figure 12, the fuel flow rate of Ul.8 seem and air flow 
rate of 2^50 seem used in the preliminary pollution studies in this report were 
selected as a reasonable compromise between a sufficiently large and yet rea- 
sonably stable flame. However, in later studies at these flow rates, an anoma- 
lous zone of instability was encountered in the pressure range from about 9 to 
16 atm. At these pressures, the flame tip pulsed or bounced vertically up and 
down but without any perceptible lateral motion. This bouncing of the flame 
tip is not believed to be caused by the controlling action of the back-pressure 
regulator because flame stability was reestablished at higher pressures; above 
15 atm, and particularly at about 20 atm, the flame was once again very stable, 
completely devoid of motion, and quite ‘majestic in its slender brilliant 
appearance. This region of stability persisted up to about 35 or Uo atm, above 
which a second type of instability sometimes occurred in the tip of the flame. 
This instability, which was more noticeable at 50 atm and is indicated in fig- 
ure 12(d), consisted of the upper portion of the flame swaying slowly from side 
to side. This swaying motion was not large, perhaps less than 1 mm, and was 
unaccompanied by vertical motion. It is possible that this instability is 
hydrodynamic, as proposed by Powell and Browne (ref. IT). At such high pres- 
sures , strong density gradients exist in the flame because of the very thin 
reaction zone and its accompanying steep concentration and temperature gradi- 
ents. Such gradients were particularly obvious just above the flame, even with 
the unaided eye, from the strong refractive index gradients appearing about one 
full flame height above the flame. Since these steep density gradients, and 
their necessarily associated steep velocity gradients, undoubtedly occur within 
the flame as well, it is not unreasonable to suspect that these gradients may 
be responsible for the observed flame instability. 


Oxides of Nitrogen 

After sufficient experience had been gained in operating the high-pressure 
flame system, experimental tests were made over the entire pressure range from 
1 to 50 atm to determine the effect of pressure on the formation of N0y. All 
these tests were made with methane as the fuel and air as the oxidizer. A 


18 



single flow condition, 4l.8 seem for fuel and 2450 seem for air,^ was studied 
for the reasons mentioned in the previous section. During the tests the stain- 
less steel filter ahead of the hack-pressure regulator was removed periodically 
and inspected for soot accumulation. None was observed because the experimental 
conditions were selected so that virtually no free soot, or smoke, was liberated 
from the flame during a test. In those few instances when the flame was inad- 
vertently ignited at too high a fuel flow rate and smoke was generated, the test 
was aborted and the filter was replaced. For all tests the filters employed 
had been precleaned ultrasonically with freon and oven-dried for 2 days at 
423 K to remove any traces of oil or contamination in the filters. 

Results of these experimental tests are shown in figure 13 as molar emis- 
sion index of NOx, I m , as a function of pressure. The molar emission index 
of NOx is defined as the number of moles of NOx (combined total concentrations 
of NO and N0 2 ) produced per mole of fuel burned. Experimentally, it was 
obtained from the ratio of the measured concentrations of NOy and C0 2 in the 
effluent from the high-pressure chamber. The molar emission index is analogous 
to the emission index I customarily used to report data on gas turbine 
engines. The emission index is defined as the number of grams of NOx produced, 
as if all were present as N0 2 , per kilogram of fuel burned. According to these 
definitions, then, the emission index for methane is 2875 times the molar emis- 
sion index. The data points shown in figure 13 for any given pressure were 
taken on separate days and indicate a remarkable degree of reproducibility. 
Figure 13 shows that the molar emission index which is 0.85 x 10-3 a t 1 atm 
increases rapidly as pressure is increased with a maximum value of 1.21 x 10 “ 3 
being reached at about 9 atm; above this pressure the molar emission index 
decreases slowly and continuously until at 50 atm it is back to roughly the 
same value that it was at 1 atm. Particularly significant for gas turbine 
engine applications is that for all pressures above 1 atm and approaching 
50 atm, the molar emission index is greater than at 1 atm. Limited tests at a 
few selected pressures in which the flow rate of methane was intentionally 
varied showed little effect on I m over the methane flow rate range from 
37-6 seem to 47.9 seem. No attempts were made to take data at subatmospheric 
pressures . 

A molar emission index for NO is not reported since little could be learned 
from such a quantity; much of the NO formed in the combustion process is oxi- 
dized to N0 2 in the postflame gases and in the sample lines leading to the 
chemiluminescent analyzer. Because total oxides of nitrogen are conserved both 
in sampling and in analysis by chemiluminescence in lean flames (ref. l8), the 
reported concentration of NOx is an accurate indication of the sum of the con- 
centrations of NO and N0 2 formed during the combustion process. In any event, 
it is the total oxides of nitrogen which are of interest, and it is these total 
oxides which are reflected by the molar emission index for NOx presented in 
figure 13. 

Literature data on NOx amissions as a function of pressure during combus- 
tion are limited, both in the number of studies and in the pressure range cov- 
ered. Nevertheless, several such studies were located, three for gas turbine 
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combustors, and one for an opposed-jet burner. It is important to recognize 
that although all these studies are for diffusion flames, they are for highly 
turbulent systems -with geometry, fuels, and inlet conditions markedly different 
from those in the present study. However, the basic trend of the data from 
these studies is of some interest in comparison with the present results. This 
comparison is shown in figure lk . Since three of these studies are for gas tur- 
bines or gas turbine combustors with the results reported as emission index (i), 
the present data were converted to emission index for direct comparison. The 
data for these gas turbine combustor studies in figure lk are plotted against 
combustor inlet air pressure over the actual pressure range of the experimental 
data. These studies are (l) Diehl and Holdeman (ref. 19) for (liquid) Jet-A 
fuel in a JT8D-9 engine combustor system with constant combustor inlet air tem- 
perature of about 659 K, (2) Norgren and Ingebo (ref. 20 ) for (gaseous) propane 
fuel in a test combustor with a constant combustor inlet air temperature of 
589 K, and (3) Quigg (ref. 2l) for a typical liquid jet fuel in an experimental 
combustor with constant inlet air temperature of about 530 K. The air tempera- 
tures in these turbine studies were considerably above the inlet air temperature 
for the present methane data, which was only about 295 K. From figure lU, the 
trend of increasing emission index with increasing inlet air temperature is 
apparent. A line of slope 1/2 has been drawn in figure lH for reference pur- 
poses, since it is commonly assumed (see, for instance, ref. 22) that the 
NOx emission index in gas turbine combustors correlates as p 1/2. 

The results of Quigg (ref. 21 ) are particularly interesting from two points 
of view. First, they are similar to the present methane data in that they 
exhibit a maximum in the emission index as a function of pressure, and secondly, 
they are representative of data for a wide variety of liquid hydrocarbons. 

Quigg conducted an extensive parametric study in a test combustor for a wide 
variety of pure liquid hydrocarbons along with JP-5 jet fuel for comparison. 
Included in his parametric study were variations in combustor inlet air tem- 
perature and humidity, combustor pressure, heat input, reference velocity, and 
method of fuel injection. Regression equations empirically relating emission 
index to the various combustor operating parameters were then formulated. For 
all fuels, with all operating parameters constant except pressure, these regres- 
sion equations are of the form 


loglO %o x = a + gp - 0.00l6p 2 


where the precise values of a and 3 are determined by the values of the 
various operating parameters and p is the combustor inlet pressure in atmo- 
spheres. Typical values for a and 3 are on the order of 0.60 and 0.0^2, 
respectively ; these values were used in generating the curve labeled !, Quigg 
(ref. 21)" in figure lk. 

The dashed curve in figure lk labeled "smoky" is from Norgren and Ingebo 
(ref. 20) for a mixture of ASTM A-l (liquid) jet fuel and (gaseous) propane in 
a different nozzle system from that employed for the propane data shown. Mix- 
ing was very poor in this nozzle system, and the- large amount of smoke produced 
indicated an overly rich local condition in the primary zone of the combustor. 
The maximum in this curve is obvious. 
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In these comparisons of turbine combustor data with the present laminar 
diffusion flame data (fig. l4), only the shape of the curve, that is, the influ- 
ence of pressure on NOy formation, is of significance. The absolute magnitude 
of the emission index depends on a large number of additional factors, for 
example, the fuels used and the operating conditions in a gas turbine combustor 
as compared with a laminar diffusion flame. However, the fact that maxima have 
also been observed in the combustion of hydrocarbons in gas turbine combustors 
as well as in the present laminar diffusion flame suggests that this phenomenon 
is worthy of further study. 

The effect of pressure on NOx formation has also been studied in flames 
other than those of hydrocarbon fuels. In an early study by Newitt and Lamont 
(ref. 23), NOx formation was measured during combustion of CO in artificial 
atmospheres of O 2 /N 2 over the pressure range from 10 to 100 atm in an opposed- 
jet diffusion flame. (The purpose of this study was to determine whether 
increased pressure could be effectively employed to produce NO efficiently. ) 

Some of the data of reference 23 have been abstracted for (J ) 1 = 0.5, and these 
are shown in figure l4 for three different supporting atmospheres of 02 /N 2 * 

The emission index for these CO flames, even In an atmosphere of composition 
fairly close to that of air (33.3 percent O 2 / 66.6 percent N 2 ), is more than an 
order of magnitude greater than the emission index measured for methane in the 
present study. This is likely due to the fact that CO flames have much higher 
adiabatic flame temperatures than CH 4 flames, about 225 K higher (ref. 24), and 
the radiation loss from CO flames should be the same or less than that from 
CHl| flames. In contrast to the hydrocarbon data previously discussed, the 
emission index for CO flames in figure l4 shows no tendency to reach a maximum, 
at least below 100 atm. Several important differences between CO flames and 
hydrocarbon flames are the higher flame temperature for CO flames and the fact 
that no hydrocarbon radicals are present and no carbon is formed in CO flames. 
Regarding this latter point, note the similarity between the gas turbine com- 
bustor data in figure l4 labeled M smoky" and the present methane diffusion 
flame data; both show an obvious peak in the emission index curve as a function 
of pressure, and both have a high degree of carbon formation. 

Certainly, any attempt to explain the maximum in the NOx molar emission 
index curves in figures 13 and l4 cannot ignore the influence of temperature 
because of the known strong temperature dependence of the NO-forming reactions. 
However, in light of the experimental temperature data presented in the follow- 
ing section, changes in flame temperature with pressure cannot be used to 
explain the shapes of these curves. 


Flame Temperature 

Estimates of the flame temperature were obtained from measurements of 
brightness temperature at 5000 A and 6500 A as discussed previously. It should 
be recognized that the flame temperatures measured by this technique are really 
average carbon temperatures , that is , spatial average temperatures of the hot 
carbon particles existing in the flame in the zone viewed, with the optical 
pyrometer(s) . It has been shown by Kuhn and Tankin (ref. 5 ) in their work with 
a large propane diffusion flame that the temperature of these carbon particles 
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is, in fact, a good indication of the true temperature of the flame gases in the 
vicinity of the carbon particles. But their results also indicate that the peak 
temperature reached in the flame may be as much as 100 K higher than that in the 
carbon zone. However, all their work was done at 1 atm where a separation of 
the carbon zone from the (blue) gaseous reaction zone must undoubtedly exist 
(they did not show color photographs of their flames such as those in fig. 7). 
Because of the much lower carbon-forming tendency of the small methane flames 
in the present study, no meaningful temperature measurements could be made at 
pressures below 5 atm; and by the time this pressure was reached, the concen- 
tration of carbon in the flame had dramatically increased, the reaction zone 
had become thinner and more condensed, and the separate blue zone was no longer 
evident. (See fig. 7*) Hence, it may be assumed that no significant separate 
gaseous reaction zone exists and that all important combustion reactions occur 
simultaneously within the visible luminous carbon zone. Since flame tempera- 
ture was measured in the tip of the flame in the region of most intense radia- 
tion, this measured temperature may be taken as a reasonable representation of 
the peak temperature occurring in the flame. 

Figure 15 shows a plot of the measured flame temperature T^ as a func- 
tion of pressure with data limits at each pressure indicated by a vertical line. 
The fairly broad limits are caused by the high sensitivity of the flame tempera- 
ture, calculated by equation (3), to small changes in the two individually 
measured brightness temperatures at the relatively close wavelengths, 5000 A 
and 6500 A. Although a wider separation between the two wavelengths would be 
desirable to minimize errors, wider separation cannot be conveniently achieved 
because the relationship a k/A a has only been shown to apply in the visible 
region. An estimate of the flame temperature at pressures below 5 atm was 
obtained by extrapolating the temperature curve to a value at 1 atm arbitrarily 
taken to be 100 K below the adiabatic flame temperature of 2230 K to account 
for molecular radiation and radiation from the small amount of carbon in the 
flame. This value of 2130 K is in agreement with the maximum value of 2116 K 
calculated from thermocouple probe data reported by Tuteja (ref. 25)^ for large, 
1-atm methane-air diffusion flames. For comparison, the adiabatic flame tem- 
perature T ad for methane-air combustion at stoichiometric conditions has also 
been plotted in figure 15 . 

Figure 15 shows that the measured temperature Tf decreases about l60 K 
from 1 atm to 5 atm, while carbon concentration, and hence radiation from the 
flame, substantially increases. There is a further decrease of about 120 K 
from 5 atm to 20 atm as carbon concentration continues to increase. (See 
figs.- 7 and 10.) From 20 atm to 50 atm, only a 50-K decrease in temperature 
takes place. There is little similarity between the shape of this temperature 


Tuteja reports a maximum thermocouple temperature of about 1978 K 
(3100° F) uncorrected for thermocouple radiation losses. Using standard radia- 
tion thermocouple corrections (ref. 26) and the usually accepted value of 0.25 
for thermocouple emissivity (see, for instance, refs. 27 and 28) result in a 
maximum corrected flame temperature of 2116 K. It may be observed that Tuteja 
himself uses an excessively large emissivity, 0.U0, which yields a maximum 
flame temperature only 20 K below the adiabatic flame temperature. 
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curve and the shape of the curve of molar emission index in figure 13- At 
pressures just above 1 atm, molar emission index increases rapidly while tem- 
perature decreases; near 9 atm, molar emission index peaks while temperature 
decreases smoothly; and above about 20 atm, molar emission index continually 
decreases while temperature decreases only slightly. Hence, it is clear that 
the shape of the temperature curve with pressure cannot explain by itself the 
shape of the curve of molar emission index with pressure. Possible explana- 
tions for the shape of this molar emission index curve are explored in some 
detail in the following section. 


INTERPRETATION OF N0 X DATA 

The fact that the molar emission index for N0 X for the present methane-air 
diffusion flame attains a maximum as a function of pressure and the fact that 
the direction of change of molar emission index with pressure is, in general, 
contrary to the direction of change of temperature with pressure are unusual 
results and warrant further consideration. In this section, experimental aver- 
age rate of formation of N0 X is estimated and compared with rates predicted for 
a variety of different conditions by the conventional two-step Zeldovich mecha- 
nism for NO formation. The experimentally observed peak in the curve of molar 
emission index as a function of pressure is discussed, and several possible 
explanations are offered. 


Rate of Formation of N0 X 

An expression for the average rate of formation of N0 X in the flame can be 
obtained from the definition of molar emission index I m . Defined as the num- 
ber of moles of N0 X produced per mole of fuel consumed, the molar emission index 
can, under conditions of steady state and complete combustion of the fuel 
(methane), be expressed as 


I 


m 


n N0 X 

n CH4 


. • 

where n^j^ is m °l ar flow rate of CH 4 into the flame and n^o^. is the 

molar flow rate of N0 X out of the flame. Rearranging this expression and 
dividing by a reaction volume yields 


n N0 X 


= I. 


n CH4 


'Rx 


m 


V- 


Rx 


If the not unreasonable assumption is made that the reaction volume in which 
N0 X is formed is approximately equal to (though not necessarily coincident with) 
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the reaction volume in which CH^ is combusted to CC >2 and H 2 0,^ this equation 
can be rewritten as 


-E -E 

r N0 x m r CH4 


(k) 


, -E 

where r. 


NCK 


is the average rate of formation of N0 X within the reaction volume 


-E 


in which it is formed and r n is the average overall rate of conversion of 


CH^ to CC >2 and H 2 O within the reaction volume in which such conversion occurs. 
The reaction rate r?L can be obtained by dividing the known molar flow rate 

CE h 

of methane n cHj^ the volume of the flame zone, which can be obtained from 

the color photographs of the flame in figure 7. Reaction volumes obtained from 
these photographs should be a fairly accurate measure of the overall fuel con- 
sumption volumes since it has been experimentally demonstrated that methane 
combustion is complete within the flame (i.e., there are no significant quan- 
tities of unburned hydrocarbons, CO, or smoke liberated from the flame). Hence, 
the average methane reaction rate so determined should be a fairly realistic 
estimate of the true average overall rate of consumption of methane in the 
flame. 


In defining an appropriate methane reaction volume V-p x from the photo- 
graphs in figure 7, both the blue and the luminous carbon zones must be 
included because no distinction should be drawn between methane consumption 
entirely in the gas phase and methane consumption when (solid) carbon is formed 
as a reaction intermediate; that is, there is no fundamental difference in the 
concept of reaction volume in which methane is consumed whether the path be 


or 


CH^ + 20 2 


Gaseous reaction ^ + 2H Q 

intermediates ^ ^ 


CHj^ + 20 2 


Gaseous reaction 
intermediates 


+ 


Solid reaction 
intermediates 


C0 2 


+ 2H 2 0 


Prior to using the photographs in figure 7 for obtaining estimates of such 
reaction volumes, however, it had to be ascertained whether all zones of 


^Certainly it is reasonable to assume that the volume in which N0 X is 
formed is roughly proportional to the volume in which CH^ is consumed. It is 
considered unlikely that the constant of proportionality would be much less 
than 1/2 or much greater than 2; hence, for purposes of convenience, it has 
been taken to be unity. In any event, the precise value selected for this 
proportionality constant within these bounds has no material effect on the 
subsequent discussion. 
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significant reaction could be seen in these photographs. In separate experi- 
ments, photographs were taken of the flame from 1 to 50 atm through interference 
filters at two wavelengths: 3100 A where OH band emission occurs and 8700 A 

where only continuum emission from carbon particles is expected. On comparing 
these photographs with the color photographs in figure 7 * it was found that the 
photographs taken with the 3100-A filter coincided with the blue and green zones 
for the lower pressure flames and with the yellow luminous carbon zone for the 
higher pressure flames, and that the photographs taken with the 8700-A filter 
coincided with the white or yellow luminous zone at all pressures. To deter- 
mine the nature of the flame emission at these separate wavelengths, an emission 
spectrum of the flame was taken with a 0 . 25 -m monochromator in the spectral 
regions of 3100 A and 8700 A at both 1 and 50 atm. At 1 atm, the OH rotational 
band structure was clearly visible in the vicinity of 3100 A, whereas at 87 OO A, 
only the continuum from carbon-particle emission was evident. At 50 atm, the 
carbon-particle continuum at 8700 A was still observed; but the OH band struc- 
ture was no longer evident near 3100 A, there being only a faint continuum of 
increasing intensity with increasing wavelength, undoubtedly due to carbon- 
particle emission. The lack of OH band structure could be due either to a very 
low OH concentration in the flame at 50 atm, or to excited OH failing to the 
ground state by collision damping rather than by emission. On the basis of 
these separate spectral studies it was concluded that the color photographs in 
figure 7 did indeed give a fair representation of the significant reaction 
zones, including both the OH region and the luminous carbon zone, and, further, 
that it was unlikely that zones of significant gas-phase reactions were missed, 
even for the higher pressure flames where carbon-particle luminosity becomes 
so intense. 


To obtain estimates of the methane reaction volumes, the reaction zones 
shown in the color photographs were traced on fine-grid graph paper, and the 
areas so defined were analytically revolved about the center line of the flame. 
In estimating such reaction volumes, possible image distortion, caused by 
refraction of light by density gradients in the flame, was ignored. A plot of 
the resulting reaction volumes as a function of pressure is shown in figure l 6 
along with their estimated error bounds. 


If the reaction volumes given in figure 1 6 are divided into the known 
(constant) flow rate for methane of 4l.8 seem expressed in moles/sec, the aver- 


age overall methane consumption rate r 


-E 


CH4 


is obtained. This rate is plotted 


as a function of pressure in figure 17. As limiting cases for methane consump- 
tion, the average overall rates of consumption of methane based on gas-phase 

— s t o ich 

" and based on carbon-particle oxidation kinetics 


reaction kinetics 
-carbon 
r CH4 

ure 17 for comparison. 


CH4 

were calculated as described subsequently and are also plotted in fig- 


tion 


In calculating the limiting reaction rate for 
r , the following assumptions are made. 


gas-phase methane oxida- 
All methane is consumed by 
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reaction with oxygen, and the concentration of oxygen in the presence of 
methane during reaction is just that needed for complete combustion, that is, 
the stoichiometric amount. For this stoichiometric combustion, the mole frac- 
tions of methane and of oxygen xq 2 are taken to be 0.095 and 0.190, 

respectively, the same amounts that would exist in a premixed stoichiometric 
methane-air flame. The kinetic expression used for calculating the overall 
rate of methane consumption (i.e., the single-step conversion of CH]| and to 
CO2 and H2O) is 


rcffl, ■ lO 1 * 1 ^OO/BTc^ 


(5) 


which is based on a study of methane-air diffusion flames by Tuteja (ref. 25). 
Because the temperature in equation (5) represents an average temperature 
throughout the reaction zone, which is somewhat uncertain but obviously smaller 
than the measured peak flame temperature T f , several average temperatures less 
than Tf were employed. These were arbitrarily taken to be 50 K, 100 K, and 

— st oi ch 

150 K below T-p, and the kinetic reaction rate r nTJ plotted in figure 17 

has been calculated for these three average temperatures. There are only minor 
differences in the predicted kinetic reaction rates for the three temperatures. 
The fact that all these kinetic rates are more than 2 orders of magnitude 

-E 

greater than the experimental average reaction rate r pw indicates that the 

experimental reaction rate is strongly diffusion controlled at all pressures, 
as expected. 


Since at high pressure, methane consumption proceeds to a large extent 


-carbon 

CH], 


has also been 


through (solid) carbon intermediates, a reaction rate r 

calculated. This reaction rate is based on the assumptions that all methane 
is instantaneously and completely pyrolyzed to carbon particles and hydrogen, 
and that oxidation of these carbon particles is the rate-controlling step. 

To estimate this oxidation rate of carbon particles, the expression of Nagle 
and Strickland-Cons table for the oxidation rate of pyrolytic graphite was 
used in accordance with the suggestion of Appleton (ref. 29) that it fairly 
accurately describes the oxidation rate of carbon particles within a flame. 
For this reaction, as before, the stoichiometric amount of oxygen was assumed 
to be in the reaction vicinity of the carbon particles (and hydrogen). Since 
the reaction-rate expression of Nagle and Strickland-Constable (as given by 
Appleton (ref. 29) and reproduced in appendix E) is per unit area of carbon 
surface, the size of the carbon particles in the flame affects the overall 
reaction rate. Calculations were made for spherical particles, both 50 1 and 

o 

100 A in diameter, sizes which are within the bounds observed by Narasimhan 
and Foster (ref. l6) in their study of carbon-particle growth in turbulent 

—carbon 

methane-air flames. For simplicity, the calculation of r OTJ was made at 

urii| 

only one average temperature, T av = Tj> - 100 K. The resulting reaction rate 

—carbon 

is plotted in figure 17. Although the kinetic prediction for r is 
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below that for » it is still more than an order of magnitude greater 

-E ^ 

than r CH U - Therefore, regardless of the reaction path for the consumption of 

methane, the reaction rate in the present diffusion flame appears to be defi- 
nitely limited by diffusion. 

The average rate of formation of nitrogen oxides r in the flame was 


NCk 


-E 


calculated from the plot of r c in figure IT, the plot of I m in fig- 

ure 13, and equation (4). A plot of r as a function of pressure is 

shown in figure 18. This estimate of the average rate of formation of NO^ 
in the flame can now be compared with formation rates predicted by the con- 
ventional two-step Zeldovich mechanism for NO formation: 


0 + n 2 J NO + N 
k r 


N + 0 2 t NO + 0 
k r 


( 6 ) 


(T) 


The generally accepted values for the reaction-rate constants, given by Bracco 
(ref. 30) , are 

k fg = 1.36 x 10 14 e - 75l * 00/RT 

k rg = 3.1 x 10 13 e - 334 /* 1 

k f7 = 6.k x 10 9 T e- fe50/RT 

k , = 1.5 X 10 9 Te- 3661,0/RT 

r7 


If elementary rate expressions, based on equations (6) and (7), are written for 
the rate of formation of NO (dCjjo/^'t)^ the usual steady-state assumption for 
nitrogen atoms is made, and an order of magnitude analysis is used to eliminate 
negligible terms, the customary expression for the rate of formation of NO is 
obtained : 


dC 


r N0 


NO 


dt 


= 2 k f6 C 0 C N. 


( 8 ) 
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One approach which has "been used in the literature in applying equa- 
tion (8) has "been to assume that oxygen atoms are in equilibrium with oxygen 
molecules : 


0 


2 



20 + M 


where the equilibrium constant K Q given by Bracco is 


K. 


^ = = o ck * m 1 * I -118000/RT 

Cq 2 “ k r9 “ ^ T 


( 9 ) 


( 10 ) 


Substitution of equation (10) into equation (8) permits the calculation of the 
predicted reaction rate, r^ Q , Provided that and are known. Unfor- 

tunately, unlike the situation for a premixed flame, there is some uncertainty 
regarding the proper values for these concentrations in a diffusion flame. 

For nitrogen, species concentration profiles in fuel-air diffusion flames 
reported by Tuteja (ref. 25) and Takagi and coworkers (ref. 28) show a large 
concentration of nitrogen, virtually the same as in air, extending well into 
the combustion zone. Hence, to a very good approximation, the mole fraction 
of nitrogen in the region of important NO formation may be taken to be 0.79. 
The proper concentration to use for oxygen is considerably more uncertain. 
Obviously, the maximum concentration of oxygen that can exist is a mole frac- 
tion of 0.21; this value, therefore, represents an upper limit. A lower limit 


for oxygen mole fraction can be obtained by replacing r 


with 


-E 

r CH, 


CH, 


in equation (5) 


from figure 17 and solving for that Cq which would give the 


observed methane reaction rate at stoichiometric concentrations of methane and 
oxygen in the flame zone. This was done for T av = T^ - 100 K, and the values 
of Co 2 obtained are plotted in figure 19 as average mole fraction xq 2 as a 

function of pressure. These values of xq^ are considered a lower limit of 

oxygen mole fraction for NO formation because they represent the average mole 
fraction of oxygen in the combustion zone, whereas maximum NO formation is 
known to occur slightly to the air side of the combustion zone where tempera- 
ture is highest and oxygen concentration is somewhat greater than in the fuel 
combustion zone. To get an intermediate, more realistic value for oxygen mole 
fraction which may actually exist in the NO formation zone of a diffusion 
flame (at least at 1 atm), species profiles reported by Tuteja (ref. 25) and 
Takagi and coworkers (ref. 28) in 1-atm diffusion flames can be used. Profiles 
reported in both references indicate that the mole fraction of NO in the flame 
has virtually the same shape as the temperature profile, which peaks slightly 
to the air side of the combustion zone. (A qualitative sketch of a typical 
profile is depicted in fig. 20.) Since NO seems to be formed in the highest 
temperature zone in the flame, a number of species profiles from references 25 
and 28 for the mole fraction of oxygen were integrated over that region of the 
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This gave an 


flame within about 300 K of the peak in the temperature profile.^ 
approximate average oxygen mole fraction of about 0.04. 

The three values for xq 2 presented in the preceding paragraph were used 
in equations (8) and (10) to predict ? N q . Such predictions are plotted in 

figure 18 for one average temperature, T = T« - 100 K. Certainly the most 

• dv i _ 

striking aspect of these predicted rates in figure 18 is that they are as much 

as 2 orders of magnitude smaller than the experimentally observed rate. To 
dispel the possible argument that T av = Tf - 100 K is really too low, similar 
rate predictions were made for XO 2 = 0.21 and 0.04 for T^, the peak tempera- 
ture observed in the flame. These rates are also plotted in figure 18, and 
they, too, fall well below the experimentally observed rate. The predicted 
rates are only weakly dependent on pressure because the decreasing flame tem- 
perature is sufficient to cancel any increase in rate caused by increasing 
pressure. The slight increase in the predicted rates between 20 atm and 50 atm 
is due largely to the increase in pressure, since flame temperatures between 
20 atm and 50 atm differ by only 50 K. In strong contrast to the shape of the 
predicted rate curves is the shape of the experimental rate curve, which 
increases markedly with pressure at low pressures but begins to level off at 
higher pressures. 

Experimental rates of NO-^ formation considerably higher than rates pre- 
dicted by the Zeldovich mechanism, assuming oxygen-atom equilibrium, have often 
been observed and have been the subject of much discussion in the literature. 
The discrepancy is generally attributed to the presence of superequilibrium 
concentrations of oxygen atoms in the flame; that is, during the combustion 
process, oxygen atoms are formed as reaction intermediates in concentrations 
greater than those which would exist if the oxygen atoms were in chemical 
equilibrium with oxygen molecules in the flame. (See, for instance, refs. 27, 
31, 32 , 33, and 34 for premixed flames and ref. 28 for diffusion flames.) 

Since the Zeldovich mechanism is still accepted as the proper reaction mecha- 
nism, at least for lean and near-stoichiometric flames, the problem then 
becomes one of experimentally determining what the true mole fraction of oxygen 
atoms really is. Customarily the mole fraction of oxygen atoms is determined 
either from a detailed kinetic combustion mechanism or by calculating that 
oxygen-atom mole fraction required to give the experimentally observed rate 
of NOx formation. For this purpose, the rate expression in equation (8) based 
on the Zeldovich mechanism is usually used,- and results are usually expressed 
as ratios of superequilibrium concentration to the equilibrium concentration 
of oxygen atoms. Such calculations have been made for the present data for 
T av = Tf - 100 K, and the results are plotted in figure 21 as the ratio 

x n , j/x n . The mole fraction x n is that oxygen-atom mole 

The value of 300 K was chosen for the following reason. If equation (10) 
is substituted into equation (8) and an average reaction temperature on the 
order of 2000 K is assumed, the rate of NO formation drops by 2 orders of 
magnitude for a 300-K temperature decrease. Hence, all significant NO forma- 
tion must occur within at least 300 K of the peak temperature. 
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fraction required "by equation (8) to give the observed experimental rate of 

stoich 


formation 


-E 

r N0 x 


for 


% = °’ T9 


and 


T = T - 100 K, and 
av f 


0 ,eq_ 


is the 


equilibrium oxygen-atom mole fraction for stoichiometric premixed combustion 
at the same temperature T av calculated with the computer program of Svehla 
and McBride (ref. 35)- 


From figure 21, the required oxygen-atom ratio is about 30 at 1 atm, 
increases to a maximum value of about 10 000 near 35 atm, and thereafter 
decreases slightly. The value of 30 at 1 atm is in the neighborhood of typical 
theoretical and experimental values reported for premixed flames at cf> = 1. 

(See table I.) Hence, the present data are not inconsistent with other results 
in the literature at 1 atm. But as pressure increases, the required oxygen- 
atom ratio increases to very large values. The only known study reporting the 
effect of pressure on required oxygen-atom ratio is that of Ay and Sichel 
(ref. 34), who theoretically calculated required oxygen-atom ratios for methane- 
air flames at <f> = 0.6. These ratios increase from about 100 at 1 atm to only 
about 180 at 30 atm. This increase by a factor of about 2 is considerably less 
than the factor of over 300 observed in the present study, but the calculations 
of Ay and Sichel (ref. 34) are for a much lower value of <f) (0.6) than that 

assumed to be effective for diffusion flames (1.0). It must also be appreci- 
ated that a diffusion flame differs considerably from a premixed flame, because 
even though it may be reasonable to assume an overall average equivalence ratio 
of unity throughout the combustion zone, in actual fact the local equivalence 
ratio decreases from very large values on the fuel side to virtually zero on 
the air side of the flame zone. (See fig. 20.) Table I shows that for iso- 
thermal combustion (refs. 31 and 33) the required oxygen-atom ratio increases 
with increasing equivalance ratio, but for adiabatic combustion (ref. 34) the 
required oxygen-atom ratio is a minimum at § = 1. The predictions of Ay and 
Sichel (ref. 34) for adiabatic combustion are believed to be more appropriate 
to the situation which exists across the flame zone in a diffusion flame, 
since for adiabatic combustion the peak temperature is reached near <j) = 1 and 
decreases for both larger and smaller values of <}> , a situation roughly approxi- 
mating that for local values of c}) throughout the flame zone in a diffusion 
flame. Thus, because for a diffusion flame the required oxygen-atom ratio is 
really an integrated average across the flame front, and because according to 
the calculations of Ay and Sichel, the ratio is a minimum for a local (j> of 
unity, the required oxygen-atom ratio would be expected to be somewhat larger 
than that for a premixed flame at (f> = 1. The present value of required oxygen- 
atom ratio of 30 at 1 atm (fig. 21), being somewhat larger than the theoreti- 
cally calculated value of Ay and Sichel (ref. 34) of 20, is at least in the 
right direction. 

At higher pressures the present values of required oxygen-atom ratio 
increase rapidly (fig. 21 ), a situation which is consistent with NO being 
formed in an increasingly fuel-rich environment as pressure increases. For 
instance, table I shows that for rich flames the required oxygen-atom ratio 
can reach values well into the thousands. This possible situation of NO being 
formed in an increasingly fuel-rich environment as the pressure increases can 
readily be rationalized in a diffusion flame. Many authors (refs. 27 , 31, 33, 
and 34) accept the concept of oxygen-atom superequilibrium as a satisfactory 
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explanation for WO formation rates in lean and near-stoichiometric flames; 
however, for rich flames a number of authors (see, for instance, refs. 33, 36, 
37, 38, and 39) reason that additional mechanisms involving nitrogen reactions 
with hydrocarbons or hydrocarbon fragments are necessary to explain observed 
NO formation rates. In a diffusion flame, there is an abundance of nitrogen 
well into the fuel side of the combustion zone to permit such reactions to 
occur (see fig. 20), and at higher pressures the rates of these reactions should 
be enhanced. Furthermore, the reaction zone thins as pressure increases, and 
the average mole fraction of oxygen molecules in the flame zone begins to 
decrease (see fig. 19) 9 and along with it, possibly, the oxygen-atom concentra- 
tion. Hence, it seems reasonable that an increasing fraction of the total NO 
formed could be formed by hydrocarbon-nitrogen reactions as pressure is 
increased, particularly at the lower temperatures which accompany the increase 
in pressure (lower temperatures favor reduced oxygen-atom concentrations). 

Takagi and coworkers (ref. 28) reached a similar conclusion in their study 
of turbulent hydrogen-air and propane-air diffusion flames at 1 atm, in which 
the temperatures were reduced by varying amounts of nitrogen dilution. Although 
they could satisfactorily explain their results for hydrogen flames by appealing 
to an oxygen-atom overshoot , for the propane-air flames they would have had to 
postulate unreasonably high oxygen-atom overshoots to account for their results 
in the reduced temperature flames. Accordingly, in agreement with a number of 
other authors, they reasoned that for premixed flames reactions involving nitro- 
gen with hydrocarbons or hydrocarbon fragments must be occurring. It may be of 
some additional interest to note that for the 1-atm turbulent propane-air dif- 
fusion flame of reference 28 issuing from a 2-mm-diameter fuel port, an overall 
average NO formation rate of 8.4 ppm/msec was obtained. This compares with an 
average rate of about 12.0 ppm/msec for the present laminar methane-air dif- 
fusion flame at 1 atm issuing from a 3 . 06-mm-diameter fuel port. (These low 
values for hydrocarbon fuels can be contrasted with the much higher value of 
67 ppm/msec reported in reference 28 for hydrogen-air flames with the same 
burner system. ) 


Discussion of Molar Emission Index Curve 

Unfortunately, none of the foregoing discussion explains the maximum in 
the molar emission index as a function of pressure shown in figure 13. Although 
several explanations for this phenomenon are possible, the present data do not 
permit a clear choice among them at this time.' One of the more likely explana- 
tions is the possibility that at higher pressures the NO concentration in the 
flame has essentially reached its equilibrium value, which decreases as pres- 
sure is increased. The computer program of Svehla and McBride (ref. 35) was 
used to calculate equilibrium mole fractions of NO for stoichiometric methane- 
air combustion at various assumed average flame temperatures; these mole frac- 
tions were then divided by the stoichiometric mole fraction of methane to 
obtain an equilibrium molar emission index (which is the same as that which 
would exist for a premixed flame). These values of I m e( ^ are plotted in 

figure 22 in comparison with the ^m data from the present study. At low 
pressures, Im,eq is far above the experimental I m and it is safe to say 
that the equilibrium NO concentration has not been reached. However, Im,eq 
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decreases rapidly with increasing pressure (recall that there is also an 
associated decrease in temperature), and at higher pressures, I m e( ^ tends to 

follow the decrease in experimentally measured The implication from this 

comparison is that at low pressures the increase in experimental I m is due to 

an actual increase in the rate of NO formation as pressure increases, whereas 
at higher pressures the observed decrease is due to the fact that the equilib- 
rium concentration of NO, which decreases with increasing pressure, has been 
reached in the flame. (The temperature effect here is small, because above 
20 atm, Tf, and therefore T av , decreases by only 50 K. ) 

Although this explanation is appealing, it cannot be accepted uncritically. 
The concept of stoichiometric combustion in the flame zone of a diffusion flame 
is only a reasonable overall approximation; and treatment of this zone as if it 
were effectively premixed, and making equilibrium calculations on this basis, 
is an obvious oversimplification of the true situation. The absolute magnitude 
of average temperatures T av has little fundamental significance because of 
the uncertain temperature distribution in the reaction zone. Hence, it would 
be inappropriate to conclude from figure 22 that the most reasonable value 
for T av is close to T^ - 50 K. However, the qualitative picture suggested 
by figure 22 is, nevertheless, somewhat compelling; at high pressures, a 
decrease in I m seems to be caused by equilibrium concentrations of NO having 
been attained. This proposition could be investigated experimentally by study- 
ing flames of CO and H2 and making similar equilibrium calculations. The tem- 
peratures of flames of these fuels should change very little with pressure 
because no solid carbon is formed to radiate heat from the flame as is the case 
with hydrocarbon flames. It would also be interesting to study NO formation 
near an electrically heated carbon surface undergoing air oxidation, and to 
vary the temperature of the surface as well as the total pressure. This would 
permit a detailed examination of temperature effects independent of pressure 
effects. Flames of premixed fuel and air could also offer useful information 
because no assumptions would have to be made regarding the equivalence ratio of 
such flames. However, premixed flames are considerably more hazardous to work 
with at high pressures than diffusion flames because of potential problems of 
flashback and detonation. 

Another phenomenon possibly occurring in the present methane-air flame 
could account for a maximum in the I m curve, but to get any information about 
it is extremely difficult. As pressure changes, the concentrations of oxygen 
molecules and atoms in the zone where most of the NO is formed may also change, 
and in an uncertain way. Figure 19 shows the average oxygen concentration in 
the flame zone required to give agreement between the experimentally observed 
methane reaction rate and an overall kinetic expression for methane oxidation 
assuming stoichiometric combustion. This figure indicates that the oxygen mole 
fraction may decrease by a factor of about 4.5 from 1 atm to 50 atm. However, 
the oxygen concentration in the zone where most of the NO is formed is uncer- 
tain. If the temperature gradient on the air side of the flame zone becomes 
steeper as pressure is increased (as a result of a thinner reaction zone), the 
high-temperature zone, necessary for substantial. NO formation, will move closer 
into the combustion zone where the concentration of oxygen is lower, and the 
rate of NO formation will decrease. However, as discussed previously, more 
NO may be formed through nitrogen-hydrocarbon reactions. Experimentally, both 
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concentration and temperature profiles are difficult to measure in the present 
system because of the poor accessibility of the flame for probing and also 
because of the difficulties in making spatially resolved measurements in such 
small flames with such thin reaction zones. Furthermore, the presence of large 
concentrations of carbon in the flame makes effective probe sampling addition- 
ally difficult. If some experimental means could be devised by which the flame 
could be experimentally probed for both composition and temperature, the results 
would be extremely valuable in interpreting the N0 X data for the present 
methane-air flame, as well as for flames of other fuels. One potentially 
attractive alternate solution for gaining information regarding composition and 
temperature profiles would be to formulate a detailed analytical model of the 
flame, accounting for both reaction kinetics and mass transport as a function 
of pressure. However, one severely limiting aspect to such an approach is the 
lack of suitable kinetic mechanisms for the formation, growth, and combustion 
of carbon particles within the flame. If such a detailed analytical model of 
the flame could be successfully formulated, however, it could provide much 
insight into many of the unexplained phenomena occurring in the present dif- 
fusion flame. 

A third factor which might be partially responsible for the maximum in the 
I m curve as a function of pressure is that significant amounts of (solid) car- 
bon are formed in the flame at increased pressures. Figure 11 shows that the 
number of carbon particles in the flame may increase by as much as 5 or 6 
orders of magnitude as pressure is increased from 1 atm to 10 atm and that 
there is still some further increase above 10 atm. The presence of these car- 
bon particles in the flame, though not altering the final composition of the 
effluent from the flame (i.e., CO 2 and HgO), indicates that the elementary 
reactions governing combustion at high pressure (when carbon is formed) are 
different from those at low pressure (when carbon is virtually absent). It 
has already been proposed that nonequilibrium concentrations of oxygen atoms 
in the flame are responsible for determining the rate of NO formation (at least 
in lean and near-stoichiometric flames). All that would be required for a 
reduction in NO formation rate would be a reduction in this nonequilibrium 
oxygen-atom concentration caused by a change in reaction mechanism. For 
instance, if fewer oxygen atoms were produced by intermediate reactions during 
the oxidation of carbon in the flame than during the homogeneous oxidation of 
gaseous methane, the formation of NO in the flame would tend to be retarded, 
and at a sufficiently high carbon concentration, a decrease in the molar emis- 
sion index would occur as carbon concentration . in the flame increased further. 
The presence of solid carbon in the flame could possibly reduce NO formation 
in another way. As an increasing number of carbon particles are formed in the 
flame at higher pressures, a rapidly increasing amount of carbon is tied up in 
these particles, with substantially less carbon present as hydrocarbon frag- 
ments and radicals. (An estimate of the number of carbon atoms in a carbon 
particle can be obtained from eq. (E3). From the two terms in the denominator 
on the right side of eq. (E3), it is calculated that only one carbon particle, 
100 & in diameter, ties up more than 50 000 carbon atoms.) Hence, if direct 
nitrogen attack on hydrocarbon radicals in the flame is an important inter- 
mediate step in the formation of NO, NO formation should be reduced when sub- 
stantial amounts of carbon are formed in the flame, simply because the con- 
centration of these hydrocarbon radicals is reduced. 
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Unfortunately, known kinetic mechanisms for the oxidation of methane (see, 
for instance, refs. 31 and 34), which predict the occurrence of superequilibrium 
concentrations of oxygen atoms in the flame, do not account for the formation, 
growth, and subsequent combustion of carbon in the flame. Hence, it does not 
appear to be currently possible to examine theoretically the effect of solid 
carbon in a methane flame on oxygen-atom superequilibrium concentration or 
hydrocarbon radical concentration. However, several experiments can be pro- 
posed which could provide some information on this effect. Some of these 
experiments are identical with those already proposed for investigating the 
possibility that a decrease in I m is caused by equilibrium concentrations 
of NO having been reached. One such experiment would be to measure the molar 
emission index for CO and H 2 flames. Since carbon is not formed in either of 
these flames, the absence or presence of a maximum in the molar emission index 
would tend to support or refute the proposal of a carbon effect. Although 
existing data already indicate that the molar emission index of CO flames does 
not attain a -maximum as a function of pressure (see data from ref. 23 in 
fig. l4), these data need verification because air was not used as the oxidizer, 
the opposed- jet burner geometry employed was markedly different from the present 
geometry, and complete burning was not always attained. As an additional, but 
experimentally difficult step, the CO and H 2 flames could be seeded with finely 
divided carbon particles to investigate their effect. Flames of hydrocarbons 
which produce carbon more readily than methane, such as ethylene, ethane, and 
propane, could also be investigated. With these flames, any effects due to car- 
bon should be observed at lower pressures, although these flames very likely 
will smoke more readily. Premixed stoichiometric flames of methane and air do 
not form carbon at 1 atm, and none is expected even at high pressures. Hence, 
such flames could also be used to investigate the effect of carbon on NO forma- 
tion. To study NO formation associated with combustion of carbon particles in 
a flame, a solid carbon body could be electrically heated in a stream of air. 
Significant NO formation would occur only in the region of high temperature 
immediately adjacent to the surface, the temperature of which could be both 
measured and controlled quite accurately. 

In the present study of NO formation, fuel and air flow rates have been 
held constant, and it has been assumed that molar emission index is relatively 
insensitive to them. With air in considerable excess, all the air that can 
possibly diffuse to the flame is available, and hence, the effect of air flow 
should be minimal. Since molar emission index is normalized by fuel consump- 
tion, its dependence on fuel flow rate should also be minimal. As long as the 
flame remains laminar, total combined flow rate is likewise expected to have 
only a small effect on molar emission index. However, in future studies, all 
these factors should be considered, and the sensitivity of molar emission 
index (and of temperature) to the various fuel, air, and combined flow rates 
should be experimentally determined. 


CONCLUDING REMARKS 

A high-pressure flame system has been designed and constructed for study- 
ing nitrogen oxide formation in fuel-air combustion. Its advantages and limi- 
tations have been demonstrated by tests with a confined laminar methane-air 
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diffusion flame over the pressure range from 1 to 50 atm. In all tests (except 
those in which stability limits were determined), methane and air flow rates 
were held constant at Ul.8 standard cm^/min (seem) and 2 U 5 O seem, respectively. 
The methane issued coaxially from a 3 . 06 -mm-diameter port into a coflowing 
stream of air contained within a 20 . 5 -nun-diameter chimney. 

A number of interesting phenomena have been observed in the behavior of 
this methane-air diffusion flame over the pressure range from 1 to 50 atm. 
Measurements of nitrogen oxide concentrations in the combustion products of 
this flame have been made, and these, too, reveal many interesting phenomena. 

On the basis of the present experimental data, analyses of these data, and 
comparisons with other experimental and theoretical results , the following 
observations and conclusions are offered: 

(1) The shape of the methane-air diffusion flame changes from a wide con- 
vex flame at 1 atm to a slender concave flame at 50 atm. The most notable 
change in shape occurs below about 5 atm. 

(2) Accompanying the change in the shape of the flames is a marked increase 
in carbon in the flame, and hence, the luminosity of the flame, as the pressure 
is increased from 1 to 50 atm. It is conjectured that this increasing amount 

of carbon in the flame is largely responsible for its change in shape. Accord- 
ingly, any mathematical description of the diffusional combustion of methane 
with air at pressures even moderately above 1 atm must properly account for the 
formation and subsequent combustion of carbon in the flame.. 

(3) Although the height of the flame changes to some extent with pressure, 
this change is not large, and with the exception of 50-atm flames, the flames 
were typically on the order of 10 mm high. 

(4) As total pressure of the flame is increased, the limits of stability 
and extinction are greatly reduced, and the limits within which the flame 
smokes are greatly increased. The most pronounced changes occur between 1 and 

5 atm. Regardless of flow rate, all flames at 50 atm were found to be luminous. 
Even at only 20 atm, a nonluminous flame could be produced only in a very lim- 
ited region of low fuel and air flow rates . 

(5) Concentrations of nitrogen oxides produced during combustion were 
measured and reported as a molar emission index, that is, the moles of nitrogen 
oxides formed per mole of methane consumed. The molar emission index, which is 
O .85 x 10-3 at 1 atm, increases rapidly as pressure is increased; it reaches a 
maximum value of 1.21 x 10“ 3 at about 9 atm; and thereafter it decreases slowly 
until at 50 atm it is back to the same value that it was at 1 atm. This maxi- 
m-urn occurs in spite of the fact that the peak temperature in the flame decreases 
from about 2130 K at 1 atm to about 1850 K at 20 atm, and thereafter, decreases 
relatively little as pressure increases. It is particularly significant that 
for gas turbine engine applications, at all pressures above 1 atm and approach- 
ing 50 atm, the molar emission index is greater than that at 1 atm. 

(6) Similar data for emission index for gas turbine engines, over a more 
limited range of pressures than that of the present study, also reach a maximum 
with pressure. 
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(7) For the methane-air diffusion flame, overall average methane reaction 
rates were estimated from the data and compared with limiting kinetic predic- 
tions. This comparison shows that the predicted reaction rates are more than 
an order of magnitude greater than the experimental rates and therefore con- 
firms that the combustion is diffusion controlled at all pressures. 

(8) Overall average nitrogen oxide formation rates were also calculated 
from the data. These rates are far greater than nitric oxide formation rates 
predicted by the conventional two-step Zeldovich mechanism with oxygen atoms 
assumed to be in equilibrium with oxygen molecules. 

(9) The concentrations of oxygen atoms required to give the observed 
nitrogen oxide formation rates were calculated according to the Zeldovich mech- 
anism and then divided by the equilibrium oxygen-atom concentration to obtain 
required superequilibrium oxygen-atom ratios. At 1 atm the value of this 
ratio, 30, is close to literature values for premixed flames, but increases 
well into the thousands above 5 atm, reaching 10 000 at 35 atm. At these higher 
pressures, it may be that an increasing amount of nitric oxide is formed 
through nitrogen-hydrocarbon reactions in addition to the Zeldovich mechanism. 

(10) The maximum observed in the curve of molar emission index as a func- 
tion of pressure could be due to several possible causes: 

(a) At low pressures the molar emission index is determined by the 
kinetic rate of formation of nitric oxide; hence, it increases with pres- 
sure because of the effect of increasing pressure on reaction rate. But 
at high pressures the equilibrium concentration of NO may be reached, and 
since this equilibrium concentration decreases with increasing pressure, 
the molar emission index also decreases. 

(b) The concentrations of oxygen molecules and oxygen atoms in the 
zone where most of the nitric oxide is formed may be changing with pres- 
sure in such a way as to produce a maximum value for molar emission index. 

(c) The increasing amount of carbon formed in the flame as pressure 
increases may alter the combustion mechanism to produce fewer oxygen 
atoms or fewer hydrocarbon radicals. At sufficiently high pressures, 
this decreasing concentration of oxygen atoms or hydrocarbon radicals 
would cause a decrease in nitric oxide formation. 

(11) The present preliminary data do not permit a definite choice among the 
above alternatives for explaining the observed maximum in molar emission index. 
Accordingly, a number of experiments which would aid in making such a choice 
have been suggested. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 1 , 1977 
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APPENDIX A 


DEPENDENCE OF PSEUDOEQUIVALENCE RATIO ON VOLUMETRIC FLOW RATE 
OF FUEL AND AIR AT STANDARD CONDITIONS 


Define 


S 


n. 


fuel 


in 


air /stoich 


(Al) 


Thus, for methane, S = 0.105- The pseudoequivalence ratio based on actual 
molar flow rates to the burner is defined as 


4> f 


l 

s 


n fuel 

A air 


x (pV/RT) 


fuel 


S (pV/RT) 


air 


( A2 ) 


But at standard conditions, 
( l)V 

pV _ K ; STP 
RT R(293) 


(A3) 


where Vg T p is the volumetric flow rate at standard conditions (l atm, 293 K) . 
Therefore , 


1 V S TP, fuel 
S V STP,air 


(AU) 
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APPENDIX B 


DILUTION CALCULATIONS 

The combustion of methane in air for <p ' less than 1 can be expressed as 

CH 1* + 0^f(°- 21 °2 + °' 79 "2) ” C °2 + 2H 2° + {fr - s)o 2 + p-(§f||)n 2 (Bl) 

and the mole fraction of water vapor produced then becomes 


x h 2 o 


2 


1 + 2 + - 2 

) + M 

0.79\ 

U' 

J ♦' V 

0 . 21 / 


(B2) 


For the usual conditions of operation in the present report ((f ) 1 = 0.162), 

XJJ 20 = 0*033, which is greater than the saturation value at 1 atm and room 

temperature of 296 K, about x h 2 0 = 0.028. At a temperature of 288 K (arbi- 
trarily about 8 K below room temperature), the partial pressure of water vapor 
is 0.0167 atm, which corresponds to a mole fraction of O.OI 67 at a total pres- 
sure of 1 atm. Thus, to prevent condensation in unheated sample lines, the 
mole fraction of water vapor must be less than or equal to O.OI 67 . This low 
mole fraction can be attained by nitrogen dilution, where 


x H 2 0 


1 + 2 + 



! ) + + n N 2 ,<iil 


(B3) 


where n ^,dil num ^ er moles of nitrogen added for dilution per mole 

of methane consumed. 

In the present study, the molar ratio of nitrogen to total reactants or 
combustion products customarily was 2.4, which, from equation (B3) with 
<j) f = 0.162, gives = 0.0098. This value corresponds to a saturation tem- 

perature of 279.9 K at 1 atm and 354 K at 50 atm. Therefore, the use of 
unheated sample lines at 1 atm and sample lines at elevated pressures heated 
to 393 K is adequate to prevent condensation. 
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APPENDIX C 


GAS CHROMATOGRAPHIC SYSTEM 

The gas chromatographic system employed in the early checkout tests 
consisted of two separate gas chromatographs in series as shown in figure 23 - 
The first chromatograph has two 3 . 658 -m x 0. 3175-cm-diameter columns of 
100 x 120 mesh Porapak Q, and the second has two 1.829-m x 0. 635-cm-diameter 
Molecular Sieve 5A columns . Carrier gas (helium) sweeps the sample through the 
sample loop and onto the Porapak Q column in the first gas chromatograph where 
the light gases (H 2 , 02 , N 2 , CO, and CH^) come through early with only minimal 
separation. - These gases are carried to the second gas chromatograph where they 
are separated and analyzed on the Molecular Sieve 5A column. Just before CO 2 , 
H 2 O and compounds with longer retention times are eluted from the Porapak Q 
column in the first chromatograph, the switching valve (SV, fig. 23) is thrown, 
and directs them to vent. This procedure permits analysis for the compounds 
with longer retention times on the Porapak Q column in the first gas chromato- 
graph but prevents contamination of the Molecular Sieve 5A column in the 
second. Pressure imbalance during switching is minimized by the regulating 
valves (RV in fig. 23) which have been preset to provide flow restrictions 
equal to those of the molecular sieve columns. The C 02 , H 2 O, and compounds 
with longer retention times are analyzed by thermal conductivity as they 
are eluted from the Porapak Q column, with the exception of hydrocarbons, which 
are analyzed with flame ionization. Fixed capillary restriction coils provide 
the proper splitting of the flows to the thermal conductivity and flame ioniza- 
tion detectors. 

Operating conditions for the two chromatographs are as follows: 



Chromatographs 


Porapak Q 

Molecular Sieve 5A 

Temperature, K, of - 



Column oven 

Programed 3, 

358 

Thermal conductivity cell 

i+T3 

383 

Injection port 

473 

398 

Transfer zone 

473 

383 

Bridge current , mA 

150 

200 

Helium carrier flow rate, seem . . . 

4l 

kl 

Flame ionization detector: 



Flame temperature, K 

473 

— 

Hydrogen fuel flow rate, seem . . . 

25 

— 

Air flow rate, seem 

500 

— 


a The oven temperature is programed to be isothermal at 353 K for 
3 min following injection, then to increase 15 K/min up to h^3 K and to 
hold at U 53 K for 6 min before recycling. 
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Switching of the effluent from the first gas chromatograph to vent to prevent 
this effluent from being directed to the second gas chromatograph is accom- 
plished at 1.75 to 1.80 min after sample injection. A 0.5-cm3 sample loop 
provides a sufficiently large sample for a good analysis without introducing 
too large a pressure pulse upon injection. An even larger sample loop could 
he used satisfactorily if desired. 



APPENDIX D 


ESTIMATE OF CARBON CONCENTRATION IN THE FLAME 

An estimate of the mass concentration of carbon in the tip of the flame 
can be obtained from the brightness temperature measurements in the tip of the 
flame and the calculated temperature of the flame. By using equation (2) in 
the text and by approximating the exponential by the first two terms of a 
Taylor series expansion, the spectral emissivity of the carbon particles in the 
flame can be approximated by 


e x * K X cL 


(Dl) 


The spectral extinction coefficient is given by Mie theory (see, for 

instance, ref. 40) as 


„ _ 36tt_ F(X) 

K X " p c X 


(D2) 


where is the density of a carbon particle, X is the wavelength, and 


F( X ) = 


n 2 x 


[n^(l - k^) + 2^ - 1 - Un^K^ 


(D3) 


and where n and k are the refractive and absorption indices in the complex 
index of refraction m = n - inK. Equations (Dl) and (D2) can be combined and 
solved for c, the mass concentration of carbon, to get 


c = 


p c Xe x 


36tt F(X) L 


(d4) 


This expression can be evaluated at X = 6 500 K (650 nm) by taking 
n = 1.57, k = 0.28 (see ref. 4l for propane soot), and = 2.0 g/cm- 
(a virtual constant for most carbon or soot particles) to obtain 


c = 3.38 x 10 


,-5 f*. 


(D5) 


To evaluate c, then, an estimate of the spectral emissivity is needed. 

This estimate can be obtained from the definition of spectral emissivity, 


E x( T f) e X - E x( T b,x) 


(d6) 
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where Tf is the flame temperature, ^ is the brightness 

wavelength A, and E^(T), the Planck function at temperature 


E X (T) 



temperature at 
T, is given by 

(DT) 


Neglecting the 1 in the denominator (valid for XT ^ 0.3) and substituting into 
equation (D 6) yields 



(D8) 


Values for c are obtained by substituting the brightness temperature data at 
6500 A and the calculated flame temperature into equation (D8) and by insert- 
ing the resulting values of into equation (D5) along with values for L 

obtained from the flame photographs in figure 7- These values of c are shown 
as a function of pressure in figure 10. Since the required temperature data 
are not available below 5 atm, the portion of the curve in figure 10 from 
5 atm down to 2.5 atm was constructed by a reasonable extrapolation of the 
brightness temperature data toward zero at 1 atm and extrapolation of the 
flame temperature data towards 2130 K at 1 atm (100 K below the adiabatic 
flame temperature of 2230 K to account for heat loss by molecular radiation). 

If spherical carbon particles are assumed, the number density of carbon 
particles Nq in the flame is related, by definition, to the mass concen- 
tration of carbon in the flame by 


c 


= N c 



PC 


Hence , 


N 


C 


6c 

"Pc^c 3 


0.95^9 

d c 3 


(D9) 


(D10) 


and Nq can be estimated from equations (D5) and (DIO) provided that infor- 
mation on the size of the carbon particles is available. 
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REACTION-RATE EXPRESSION FOR OXIDATION OF SOLID CARBON 

According to reference 29 , a suitable expression for the specific oxida- 
tion rate of solid carbon within a flame is given by the equation (attributed 
to Nagle and Strickland-Constable ) : 


w 


' k AP x 0 2 \ 

p * Z P*oJ + ^<1 - X) 


where w 


is in moles /cm^-sec 


and 


X = 



1 

km 

1 + P ^ 

k BP x 0 2 

2Qe -30000/RT 

4.46 x 10-3 e -15200/RT 

1.51 X lo 5 e - 9T000/RT 
21>3e Ul00/RT 


(El) 


To express the specific surface reaction rate w given in equation (El) 
as a reaction rate per unit volume in the flame, w must be multiplied by a 
factor expressing the total surface area of the carbon particles in the flame 
per unit volume of the reaction zone in the flame. For spherical particles, 
this factor is 


r 


Trd c 



(E2) 


where d^ is the particle diameter and Nq is the number density of particles 
in the flame. For the limiting case when all the carbon atoms in the methane 
fuel appear in the flame as carbon particles. 


(N c ) 

V '-/max 



(E3) 
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-where (Nq) is the number density of carbon particles in the flame zone if 

all carbon appears as soot particles; the first term in the denominator on the 
right side is the volume of a carbon particle of diameter d c ; the second term 
in the denominator is the number of carbon atoms per unit volume of a carbon 
particle, where p q is the density of a carbon particle (2,00 g/cm^), Mq is 
the molecular weight of carbon, and Na is Avogadro T s number; and N° is 

the number of carbon atoms available in the methane fuel per unit volume of the 
reaction zone in the flame assuming stoichiometric combustion in the flame 
zone. Expressing N° in terms of mole fraction yields 



x CH 1| pN A 

RT 


(E4) 


where x2tj is the mole fraction of methane fuel that would exist in the flame 

zone under stoichiometric conditions if none had been converted to soot. Sub- 
stituting equations (E3) and (E4) into (E2) and taking Xqjj = 0.095 results 
in ^ 


r max * a c T 


(E5) 


This expression is the maximum possible total surface area of carbon particles 
of diameter in the flame per unit volume of the reaction zone in the 

flame. 
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TABLE I.- REQUIRED OXYGEN-ATOM RATIOS FOR PREMIXED FLAMES 


Author 

Method 

Theoretical 

calculation 

System 

ch 4 -o 2 -n 2 

P» 

atm 

T, 

K 


x 0 ,req f d/ x 0 , eq 

Bowman (ref. 31 ) 

i 

2200 

1 

0.8 

1.0 

1.4 

o o o 

00 GO LT\ 
OJ 

Ay and Sichel 

Theoretical 

CH. -0 -N 

1 

Adiabatic 

0.6 

100 

(ref. 34) a 

calculation 

4 d d 

1 



.8 

35 







1.0 

20 




? 



1.4 

1000 




30 



. 6 

180 

Thompson ? Brown , 

Experimental 

CH, -air 

1 

1800 to 2000 

0.75 

12 

and Beer 


4 

1 

1800 to 2000 

.96 

4 

(ref. 32) 








Iverach > Basden , 

Experimental 

Ethylene- 

1 

2045 

0.8 

7 

and Kirov 


air 

1 



1.0 

20 

(ref. 33) 



I 



1.4 

1100 




▼ 

\ 


1.55 

5500 


Values of required oxygen-atom ratio attributed to Ay and Sichel were 
obtained from cross plots of figures published by them. 
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Figure 1.- Burner assembly and detail. All dimensions are in millimeters. 
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Figure 2.- Types of Burke and. Schumann diffusion flames. 
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Figure 3.- Schematic diagram of "burner and sample collection system. 
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Figure t.- Hot-wire igniter probe. Dimensions are in millimeters. 
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Figure J.~ Photographs of methane-air diffusion flames for various ambient pressures. 
Methane flow rate, 4l.8 seem; air flow rate, 24^0 seem; burner diameter, 3.0 6 mm. 
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Figure 8.- Flame height as function of pressure. 
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Figure 12.- Stability limits. 
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Figure 12.- Continued. 
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Figure 12.- Continued. 
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Figure 12.- Concluded. 
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Figure l4.» Emission index as function of pressure. 
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Figure IT.- Average methane reaction rate as function of pressure 
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FC - Flow controller 
GSV - Gas sample valve 

RV - Regulating valve serving as adjustable restriction 

SV - Switching valve 

TC - Thermal conductivity cell 

FID - Flame ionization detector 

FCR - Fixed capillary restriction coil 

Figure 23.- Gas chromatograph analytical train. 
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